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High- Tc cuprates differ from conventional superconductors in three crucial aspects: the super- 
conducting state descends from a strongly correlated Mott-Hubbard insulator (as opposed to a 
Fermi liquid), the order parameter exhibits d-wave symmetry and fluctuations play an all impor- 
tant role. We formulate an effective theory of underdoped cuprates within the pseudogap state by 
taking the advantage of these unusual features. In particular, we introduce a concept of "pairing 
protectorate" and we seek to describe various phases within this protectorate by phase disordering 
a d-wave superconductor. The elementary excitations of the protectorate are the Bogoliubov-de 
Gennes quasiparticles and topological defects in the phase of the pairing field - vortices and an- 
tivortices - which appear as quantum and thermal fluctuations. The effective low energy theory 
of these elementary excitations is shown to be, apart from intrinsic anisotropy, equivalent to the 
quantum electrodynamics in (2-1-1) space-time dimensions (QED3). A detailed derivation of this 
QED3 theory is given and some of its main physical consequences are inferred for the pseudogap 
state. As the superconducting order is destroyed by underdoping two possible outcomes emerge: 
(i) the system can go into a symmetric normal state characterized as an "algebraic Fermi liquid" 
(AFL) before developing antiferromagnetic (AF) order, or (ii) a direct transition into the AF state 
can occur. In both cases the AF order arises spontaneously, through an intrinsic "chiral" insta- 
bility of QED3/AFL. Here we focus on the properties of AFL and propose that inside the pairing 
protectorate it assumes the role reminiscent of that played by the Fermi liquid theory in conven- 
tional metals. We construct a gauge invariant electron propagator of AFL and show that within the 
expansion it has a non- Fermi liquid, Luttinger-like form with positive anomalous dimension 
r)' = 16/37r^A'^, where A'^ denotes the number of pairs of nodes. We investigate the effects of Dirac 
anisotropy by the perturbative renormalization group analysis and find that the theory flows into 
isotropic fixed point. We therefore conclude that, at long lengthscales, the AFL is stable against 
anisotropy. 



I. INTRODUCTION 

In the classic tradition of condensed matter physics, 
the phenomenon of superconductivity is usually de- 
scribed as an instability of a normal metal towards 
a quantum state in which electrons bind into Cooper 
pairs jl]. This traditional picture, further fortified 
within the Eliashberg formalism j2j, has met with spec- 
tacular success when applied to conventional, low T^- 
superconductors. This success is no accident: the BCS- 
Eliashberg theory uses the Landau theory of Fermi liq- 
uids as the underlying description of a normal metal. In 
turn, the Landau theory of Fermi liquids is one of the 
most successful theories in physics - by starting from the 
free gas of fermions it exploits the constraint on the phase 
space imposed by the Pauli exclusion principle to provide 
a comprehensive understanding of a low energy behav- 
ior of an interacting system. Among its many pleasing 
features, the Landau theory allows for methodical under- 
standing of the seeds of its own destruction - while Fermi 
liquid is hardly the true ground state of any system, |^ 
we routinely seek to understand such ground states by in- 
vestigating various instabilities of the Fermi liquid (FL) 
description, either in the particle-hole or particle-particle 
channel, few notable examples being spin- and charge- 



density waves, itinerant ferromagnetism and, of course, 
superconductivity. This remarkable success story of the 
BCS-Eliashberg-FL theory gave birth to the traditional 
paradigm: "one must understand the normal state before 
one can understand the superconductor" ||. 

The discovery of high- Tc superconductors (HTS) and 
the subsequent efforts to decipher their mysteries altered 
the above state of affairs dramatically. As experimen- 
tal information began to pour in, it became rapidly clear 
that HTS depart qualitatively from the BCS-Eliashberg- 
FL orthodoxy in at least three important ways. First, 
the high-Tc superconducting copper oxides are strongly 
interacting systems in which correlations play an essen- 
tial role. Their parent compounds, believed to contain 
one hole per copper within Cu02 layers, are far from 
good metals that the simple band-structure theory would 
predict and are instead insulating Neel antiferromagnets. 
Even when doped and the long range antiferromagnetic 
order had subsided, these materials do not seem to fol- 
low the dictates of the Fermi liquid theory, except per- 
haps in the heavily overdoped regime. Rather, the "nor- 
mal" metal state of the cuprates is anything but, and has 
been routinely dubbed "anomalous" or "strange" . In line 
with the traditional paradigm, a major theoretical effort, 
mainly inspired by Anderson [H , has been directed at first 
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understanding this anomalous "normal" state as arising 
from the Hubbard-Mott-Neel antiferromagnetic insulator 
at half-filling once a small density of mobile dopant holes 
has been induced in the Cu02 layers. The physics of 
such "doped Mott insulator" and particularly the micro- 
scopic mechanism through which the high temperature 
superconductivity itself is generated from within such a 
"normal" state, even after years of concentrated theoret- 
ical and experimental onslaught, remain a deep challenge 
to the condensed matter physics community. 

The second departure from the BCS-Eliashberg-FL or- 
thodoxy has by now been firmly established experi- 
mentally: the superconducting state in cuprates is it- 
self of unconventional symmetry Instead of a usual 
"s-wave" gap function whose magnitude generally varies 
over the Fermi surface but whose sign does not, the 
cuprates are "d-wave" superconductors. It is believed at 
the present time that this d-wave symmetry reflects the 
presence of correlations caused by strong on-site repul- 
sion and is thus a close dynamical relative of the antifer- 
romagnetic state at half-filling. Majority of authors have 
focused on such purely electronic "mechanism" of su- 
perconductivity arising from strong correlations in sharp 
distinction to the phonon-mcdiatcd pairing of the tra- 
ditional BCS-Eliashbcrg approach. Such unconventional 
symmetry of the superconducting state and the ensuing 
presence of low-energy fermionic excitations near the gap 
nodes result in a rich phenomenology of cuprates. Most 
remarkably, however, this phenomenology, at low ener- 
gies and temperatures and deep inside the superconduct- 
ing state, seems to fit within the theoretical mold of a 
model BCS-like d-wave superconductor, quite separately 
from the detailed microscopic mechanism that is at its 
origin. We consider this an experimental fact of crucial 
significance which for the rest of this paper we intend to 
fully exploit to our advantage. 

Finally, more than the intrinsic gap symmetry dis- 
tinguishes HTS from their conventional kin. They are 
also strongly fluctuating systems While strong fluc- 
tuations away from a simple mean-field description are 
familiar occurrence in magnets or liquid crystals they 
are a relatively novel phenomenon in superconductors. 
In most conventional superconductors, well-described by 
the BCS-Eliashberg-FL theory, the fluctuations are sim- 
ply not an issue: the dimensionless parameter which con- 
trols the deviations from the mean-field theory, the in- 
verse of the product of the BCS coherence length and 
the Fermi wavevector kp, is typically as small as 10"'^ 
or 10~^ and the fluctuation effects are rarely observable 
1^. In contrast, early experiments on HTS clearly es- 
tablished strong fluctuations in numerous physical quan- 
tities and various estimates put (l/S^okp) ~ 10~^ or 
larger. The strongly fluctuating nature of the supercon- 
ducting state is particularly pronounced in underdoped 
cuprates and is rather vividly manifested in recent exper- 
iments on Nernst effect by the Princeton group . The 
Nernst effect is routinely used to detect superconducting 
vortex fluctuations and a strong signal observed up to 



^Nernst ^ ^c, with TNernst Comparable to the pseudogap 
temperature T*, is most naturally interpreted in those 
terms. Other experiments have also provided both direct 
|lT| and indirect [|2j support for the pairing fluctuations 
far above Tc- The experimental evidence for such fluctu- 
ations has been compiled in Ref. fl^ . 

Recently, a different path toward the theory of HTS 
was proposed in Refs. A key aspect of this "in- 

verted" approach is the observation that, when compared 
to its neighbors in the phase diagram of cuprates, the su- 
perconducting state appears to be the "least correlated" 
and its fermionic excitations best defined A tanta- 
lizing theoretical feature of this approach is that it holds 
promise to turn the above triple predicament of strong 
correlations, fluctuations and unconventional symmetry 
in HTS into an advantage by using the superconducting 
state as its departure point. The first step in constructing 
a theory based on this "inverted" philosophy is to assume 
that the most important effect of strong correlations at 
the basic microscopic level is to build a large pseudogap 
of d-wave symmetry, which is predominantly pairing in 
origin, i.e. arises form the particle-particle (p-p) chan- 
nel. In this regard, we are following the wisdom of the 
FL theory but are replacing the free electron gas starting 
point with the free Bogoliubov-deGennes (BdG) d-wave 
quasiparticles. Under the umbrella of this large d-wave 
pseudogap, the low energy fermionic excitations enjoy 
remarkably sheltered existence. They are completely im- 
pervious to weak residual short range interactions left 
over after the effect of the pseudogap had been built in. 

There is an important new element in this parallel 
with the FL. Our reference state being a superconductor, 
we must also consider interactions of BdG quasiparticles 
with relevant collective modes of the pairing pseudogap, 
i.e. fluctuating thermal and quantum vortex-antivortex 
pairs. Within the theory of Ref. this interaction 

is represented by two U(l) gauge fields, and Ou. 
describes Doppler shift in quasiparticle energies and 
has been studied in Refs. ||l7, Its effect on low 

energy fermions is rather modest since gains mass 
from fermions both in a superconductor and in a phase- 
incoherent pseudogap state. In contrast, the Berry gauge 
field , minimally coupled to fermions and encoding the 
topological frustration inflicted upon BdG quasiparticles 
by fluctuating vortices, is massless in the pseudogap state 
and is the main source of strong scattering at low ener- 
gies . The effective theory was found to take the form 
equivalent to (2-f l)-dimensional (anisotropic) quantum 
electrodynamics (QED3) In its symmetric phase, 

QED3 is governed by the interacting critical point leading 
to a non-Fermi liquid behavior for its fermionic excita- 
tions. This "algebraic" Fermi liquid (AFL) displaces 
conventional FL as the underlying theory of the pseudo- 
gap state. 

The AFL (symmetric QED3) suffers an intrinsic insta- 
bility when vortex-antivortex fluctuations and residual 
interactions become too strong. The topological frustra- 
tion is relieved by the spontaneous generation of mass 
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for fermions, while the Berry gauge field remains mass- 
less. In the field theory literature on QED3 this instabil- 
ity is known as the dynamical chiral symmetry breaking 
(CSB) and is a well-studied and established phenomenon 
[ fi9[ , although few clouds of uncertainty still hover over 
its more quantitative aspects In cuprates, the re- 

gion of such strong vortex fluctuations corresponds to 
heavily underdoped samples and the CSB leads to spon- 
taneous creation of a whole plethora of nearly degener- 
ate ordered and gapped states from within the AFL . 
An important check on the internal consistency of the 
"inverted" approach is that the manifold of CSB states 
contains an incommensurate antiferromagnetic insulator 
(spin-density wave (SDW)) Remarkably, both 

the "algebraic" Fermi liquid and the SDW and other CSB 
insulating states arise from the one and the same QED3 
theory ||l5t] , echoing the satisfying features of Fermi liquid 
theory in conventional metals. It therefore appears that 
the "inverted" approach can be used to advance along 
the doping axis of the HTS phase diagram (Fig. 1) in 
the "opposite" direction, from a d-wave superconductor 
all the way to an antiferromagnetic insulator at very low 
doping, the low energy physics held under overall control 
of the symmetric QED3 (AFL). 



In this paper we first present a detailed derivation of 
the QED3 theory of the pairing pseudogap state in un- 
derdoped cuprates previously introduced in Ref. and 
then embark on a systematic exploration of its fermionic 
excitation spectrum and other related properties. To 
keep the paper at manageable length we confine our- 
selves to the the chirally symmetric phase of QED3, i.e. 
to the algebraic Fermi liquid (AFL) and its main proper- 
ties. The chiral symmetry broken phase is discussed sep- 
arately, in Part II of this paper. Section II and Appendix 
A contain a step-by-step manual on vortex-quasiparticle 
interactions and how the low-energy physics of such in- 
teractions can be given its mathematical formulation in 
the language of the QED3 effective theory. In Sections III 
and IV we focus on the AFL, give a detailed accounting of 
vortex-quasiparticle interactions within QED3 and com- 
pute quasiparticle spectral properties within the pseudo- 
gap symmetric phase. In Section V we then discuss the 
effects of Dirac anisotropy on the QED3 infrared fixed 
point and demonstrate that, to leading order in 
expansion, the anisotropic QED3 scales to an isotropic 
limit. Finally, we present a brief summary of our results 
and conclusions in Section VI. 



In Part II of the paper, to appear separately, we intro- 
duce the concept of chiral symmetry within the d-wave 
pairing pseudogap state pT 



22 1 . We enumerate differ- 



ent physical states within the chiral manifold and discuss 
in depth various patterns of chiral symmetry breaking 
(CSB) (Fig. 1) and fermion mass generation within our 
QED3 theory. This is followed by the applications to the 
phase diagram of cuprates within the pseudogap state. 



II. PRELIMINARIES AND 
VORTEX-QUASIPARTICLE INTERACTIONS 

A. Protectorate of the pairing pseudogap 

Our starting point is the assumption that the pseudo- 
gap is predominantly particle-particle or pairing (p-p) in 
origin and that it has a (1^2 _y2 symmetry. This assump- 
tion is given mathematical expression in the partition 
function: 

Z = j V-^\v,t) j V^{r,T) j X'(p(r,T) expos'], 

S = j dr j + + (1/5)A*A}, (1) 

where r is the imaginary time, r = {x,y), g is an ef- 
fective coupling constant in the Ax2-y2 channel, and 
= (i/j|,^/'|) are the standard Grassmann variables. 
The effective Hamiltonian 7i is given by: 



He 

A* 



A 



(2) 



with He = — f A)^ — Ei?, P = — «V (we take h = 1), 

and A the d-wave pairing operator [ESl M , 



(3) 



1 i 
^{Px, {Py, A}} - -^A{d^dy(p) 

F F 



where A(r, r) = | A| exp(i(/3(r, r)) is the center-of-mass 
gap function and {a, 6} = {ab + ha)/2. As discussed in 
Ref. the second term in Eq. (||) is necessary to pre- 
serve the overall gauge invariance. Notice that we have 
rotated A from dx2-.y2 to Axy to simplify the continuum 
limit. J 'Dlp{y^ t) denotes the integral over smooth ("spin 
wave") and singular (vortex) phase fluctuations. Ampli- 
tude fluctuations of A are suppressed at or just below 
T* and the amplitude itself is frozen at 2A ~ 3.56r* for 
T < T*. 

The fermion fields ip^ and tpi appearing in Eqs. (0J^) 
do not necessarily refer to the bare electrons. Rather, 
they represent some effective low-energy fermions of the 
theory, already fully renormalized by high-energy inter- 
actions, expected to be strong in cuprates due to Mott- 
Hubbard correlations near half- filling Q. The precise 
structure of such fermionic effective fields follows from a 
more microscopic theory and is not of our immediate con- 
cern here - we are only relying on the absence of true spin- 
charge separation which allows us to write the effective 
pairing term (^ in the BCS-like form. The experimen- 
tal evidence that supports this reasoning, at least within 
the superconducting state and its immediate vicinity, is 
rather over whelming ||, 0, Further more. Tires 

represents the "residual" interactions, i.e. the part dom- 
inated by the effective interactions in the p-h channel. 
Our main assumption is equivalent to stating that such 
interactions are in a certain sense "weak" and less impor- 
tant part of the effective Hamiltonian H than the large 
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FIG. 1: A schematic representation of tlie phase diagram of a 
cuprate superconductor. Below T* (~ TNornst), the symmet- 
ric QED3/AFL replaces Fermi liquid as the effective low en- 
ergy theory. CSB denotes chiral symmetry broken states, the 
most prominent among which is an incommensurate antifer- 
romagnet (SDW). Panel (a) represents the path between dSC 
and an antiferromagnet via the intervening algebraic Fermi 
liquid (AFL) ground state. Panel (b) shows a direct dSC- 
CSB transition with AFL describing T <C T* behavior. 

pairing interactions already incorporated through A. As 
we progress, this notion of "weakness" of Hres will be de- 
fined more precisely and with it the region of validity of 
our theory. 

The above discussion reveals that our theory is at least 
partly of phenomenological character and it must rely on 
a more microscopic description to fully define its basic 
form specified by Eqs. While it is still far from 

established just what such more microscopic description 
might be in the case of cuprates, various gauge theories 
of the t — J and related models , could all be used for 
this purpose at the present time. The main role played 
by such theories is providing reliable values of param- 
eters that feed into the basic formulation (|l|,||). These 
include, most importantly, the pairing pseudogap A it- 
self, the microscopic values of vortex core energies, the 
strengths of residual interactions, etc. Once these param- 
eters have been supplied through such an external input 
it is our task to solve for the low energy A), long dis- 
tance physics of (|l|,||). At present, these needed parame- 
ters cannot be computed reliably from a fully microscopic 
approach. We therefore combine the available theoreti- 
cal arguments with the experimentally determined phase 
diagram (Fig. 1) and argue that the pseudogap is indeed 
pairing in origin and that the transition from the super- 
conductor to the pseudogap state must proceed via ther- 
mal and quantum unbinding of vortex-antivortex pairs. 



Our "inverted" approach is similar in spirit to the Lan- 
dau theory of Fermi liquids as applied to conventional 
metals. In Fermi liquid (FL) theory the reference state 
is that of a non-interacting gas of fermions. As the inter- 
actions are turned on adiabatically, the Pauli principle 
severely restricts the available phase space for scattering 
and many of the general features of free fermion system 
are preserved, albeit in a renormalized form Q. In our 
case, the reference state is a non-interacting system of 
Bogoliubov-deGennes (BdG) fermions. One can think of 
the pairing pseudogap A as being our "Fermi energy" 
and the highest energy scale in the problem. While our 
theory cannot account for the physics at energies (or tem- 
peratures) higher than A, we will endevour to show that 
the low energy physics can be computed systematically 
and is parametrized by a handful of material constants 
whose values can be extracted either from experiments 
or from a more microscopic theory. 

In this context, the residual interactions among the 
BdG fermions, both those arising from the p-p channel 
through weak amplitude fluctuations of A and the p-h 
interactions generated by the effective spin fluctuations 
within the pseudogap state |Q, can be thought of as 
our version of "Landau interaction parameters" {-F2£+i}- 
These interactions are generically short-ranged and are 
even less effective in disturbing the coherence of nodal 
BdG fermions than their FL counterpart. With all of the 
Fermi surface gapped apart from the nodal regions, such 
interactions are irrelevant in the renormalization group 
(RG) sense and can be absorbed into renormalizations of 
various relevant parameters (more on this in Section III) . 
As the system makes the transition from a phase-coherent 
superconductor to a phase-incoherent pseudogap state, 
the nodal BdG fermions interacting only through such 
residual interactions would hardly notice - the nodal 
BdG liquid is a better Fermi liquid than the original. 

There is, however, another source of interactions 
among BdG quasiparticles which ruins this happy state 
of affairs. Our reference state being a superconductor, 
we must consider the interactions of BdG fermions with 
collective modes of A relevant near the superconductor- 
"normal" boundary, i.e. fluctuating vortices. As we 
demonstrate next, these interactions between BdG quasi- 
particles and fluctuating vortex-antivortex excitations 
play the central dynamical role in our theory. 



B. Phase fluctuations, vortices, quasiparticles and 
topological frustration 

The global U(l) gauge invariance mandates that the 
partition function (p must be independent of the overall 
choice of phase for A. We should therefore aim to elimi- 
nate the phase (/^(r, t) from the pairing term (^) in favor 
of d^if terms [/i = {x,y,T)] in the fermionic action. For 
the regular ( "spin- wave" ) piece of </? this is easily accom- 
plished by absorbing a phase factor exp(ii(^) into both 
spin- up and spin-down fermionic fields. This amounts to 
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"screening" the phase of A(r,T) (or an "XY phase", as 
commonly known) by a "half-phase" field (or "half-XY 
phase") attached to tp^ and tpi. However, as discussed 
in Ref. j2^, |2^, when dealing with the singular part of 
If, such transformation "screens" physical singly quan- 
tized hc/2e superconducting vortices with "half-vortices" 
in the fermionic fields. Consequently, this "half-angle" 
gauge transformation must be accompanied by branch 
cuts in the fermionic fields which originate and termi- 
nate at vortex positions and across which the quasiparti- 
cle wavefunction must switch its sign. These branch cuts 
are mathematical manifestation of a fundamental physi- 
cal effect: in presence of vortices, which are topological 
defects in the phase of the Cooper pair field and thus nat- 
urally bind the elementary flux of hc/2e, the motion of 
quasiparticles is topologically frustrated, since their natu- 
ral elementary flux is twice as large (hc/e). The physics 
of this topological frustration is at the origin of all non- 
trivial dynamics discussed in this paper. 

Dealing with branch cuts in a fluctuating vortex prob- 
lem is a rather cumbersome affair due to their non-local 
character and defeats the original purpose of reducing the 
problem to that of fermions interacting with local fluctu- 
ating superflow field, i.e. with 9^</5. Instead, in order to 
avoid the branch cuts, non-locality and non-single valued 
wavefunctions, we employ the singular gauge transforma- 
tion devised in Ref. |23], hereafter referred to as 'FT' 
transformation: 

exp (-i(^A) V'T 1 V'i exp {-iipB)^-'i , (4) 

where ipA + (ps — 'P- Here (pA{B) is the singular part of 
the phase due to A{B) vortex defects: 



where the transformed Hamiltonian Ti.' is: 



V X VipAiB) = 27rz ^ QiS (i 



MB), \ 



(5) 



with Qi = ±1 denoting the topological charge of i-th vor- 
tex and r^^^-*(r) its position. The labels A and B rep- 
resent some convenient but otherwise arbitrary division 
of vortex defects (loops or lines in (p{r, r)) into two sets, 
although we will soon discuss many virtues of the sym- 
metrized transformation (^ ) which apportions vortex de- 
fects equally between sets A and B . The transfor- 
mation (|^) "screens" the original superconducting phase 
(fi (or "XY phase") with two ordinary "XY phases" (pA 
and ifB attached to fermions. Both ipA and ifB are them- 
selves perfectly legitimate phase configurations of A(r, t) 
but simply with fewer vortex defects. The key feature of 
the transformation (^) is that it accomplishes "screen- 
ing" of the physical hc/2e vortices by using only "whole" 
(i.e., not "halved") vortices in fermionic fields and thus 
guarantees that the quasiparticle wavefunctions remain 
single-valued. The topological frustration still remains, 
being the genuine physical effect of the branch cuts, but 
is now incorporated directly into the fermionic part of 
the action: 

C = ij^[dr + i{drPAm + ^l[9r + i{dr(pB)]^l + 

(6) 



D 



D 



^(^-V)2 + 



2m 



with D — (Ao/2fc|,)(-n'2,7ry + TTyiTx) and t: 
The singular gauge transformation ([ 
"Berry" gauge potential 

a^, = \{d^p:>A- df^ipB) , 



p -I- a. 

generates a 



(7) 



which describes half-flux Aharonov-Bohm scattering of 
quasiparticles on vortices, couples to BdG fermions 
minimally and mimics the effect of branch cuts in 
quasiparticle- vortex dynamics [|3[ |2j, Closed 
fermion loops in the Feynman path-integral representa- 
tion of (^ acquire the (—1) phase factors due to this 
half-flux Aharonov-Bohm effect just as they would from 
a branch cut attached to a vortex defect. The topologi- 
cal frustration is now implemented through the fact that 
dfj_ipA(B) 1 minimally coupled to a loop of BdG fermions in 
the partition function (|l|) as it winds around the imag- 
inary time direction, generates the phase equal to the 
advance in the XY phase (^^(B)(r,T), which inhabits the 
same spacetime, along the closed path coinciding with 
the said fermion loop. Naturally, this advance must be 
an integer multiple of 2n - the overall factor of ^ in the 
definition of (^ reduces this further to an integer mul- 
tiple of TT. In the end, the phase factors of such fermion 
loops, which are the gauge invariant quantities, are all 
equal to ±1 and are furthermore precisely what we would 
have obtained had we opted for the "half-angle" gauge 
transformation and implemented the branch cuts. 

The above (±1) prefactors of the BdG fermion loops 
come on top of general and everpresent U(l) phase fac- 
tors exp(i(5), where the phase S depends on spacetime 
configuration of vortices. These U(l) phase factors are 
supplied by the "Doppler" gauge field 



-idf,ipA + df,ipB) 



(8) 



which denotes the classical part of the quasiparticle- 
vortex interaction. The coupling of to fermions is the 
same as that of the usual electromagnetic gauge field Afj, 
and is therefore non-minimal, due to the pairing term in 
the original Hamiltonian Ti. (|^). It is this non- minimal 
interaction with w^, which we call the Meissner coupling, 
that is responsible for the U(l) phase factors exp{i6). 
These U(l) phase factors are "random", in the sense that 
they are not topological in nature - their values depend 
on detailed distribution of superfluid fields of all vortices 
and "spin-waves" as well as on the internal structure of 
BdG fermion loops, i.e. what is the sequence of spin-up 
and spin-down portions along such loops. In this respect, 
while its minimal coupling to BdG fermions means that, 
within the lattice d-wave superconductor model [ p4| , one 
is naturally tempted to represent the Berry gauge field 
as a compact U(l) gauge field, the Doppler gauge field 
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is decidedly non-compact, lattice or no lattice p7[ . 
The reader should be advised, however, that the issue 
of "compactness" of versus "non-compactness" of 
constitutes a moot point: and as defined by Eqs. 

are not independent since the discrete spacetime 
configurations of vortex defects {vi{T)} serve as sources 
for both. We will belabor this important issue in the next 
subsection. 

For now, note that all choices of the sets A and B in 
transformation (^) are completely equivalent ~ different 
choices represent different singular gauges and w^, and 
therefore exp(j(5), are invariant under such transforma- 
tions, a^^, on the other hand, changes but only through 
the introduction of (±) unit Aharonov-Bohm fluxes at 
locations of those vortex defects involved in the transfor- 
mation. Consequently, the Z2 style (±1) phase factors as- 
sociated with that multiply the fermion loops remain 
unchanged. We now symmetrize the partition function 
with respect to this singular gauge by defining a gen- 
eralized transformation fl) as the sum over all possible 
choices of A and B, i.e., over the entire family of sin- 
gular gauge transformations. This is an Ising sum with 
2^1 niembers, where Ni is the total number of vortex 
defects in 93(r,T), and is itself yet another choice of the 
singular gauge. The many benefits of such symmetrized 
gauge will be discussed shortly but we stress here that 
its ultimate function is calculational convenience. What 
actually matters for the physics is that the original Lp be 
split into two XY phases so that the vortex defects of 
every distinct topological class are apportioned equally 
between ipA and ips ^ [§3[ |§- 

The physics behind this last requirement can be in- 
tuitively appreciated as follows: imagine that we simply 
prohibit singly quantized vortex excitations and replace 
them in by doubly quantized (hc/e) vortices. In this 
case, the Berry gauge field attaches a multiple of a full 
Aharonov-Bohm flux to each vortex position and all the 
phase factors in front of fermion loops equal unity ~ this 
is equivalent to = and complete absence of topo- 
logical frustration. It is then natural to select the gauge 
which eliminates altogether. This is straightforwardly 
accomplished by screening each doubly quantized {hc/e) 
vortex defect with a unit-flux Dirac string in -0| and 
an equivalent one in The resulting sets of A and 
B vortices are two identical replicas of each other and 
fl/i = \{^^J.'^PA — d^ifs) = 0. Now, while still staying 
within the same gauge, we allow doubly quantized vor- 
tices to relax into energetically more favorable configura- 
tions - they will immediately decay into singly quantized 
vortices and our sets A and B will end up containing an 
equal number of singly quantized vortices of each distinct 
topological class. For example, in the case of 2D thermal 
fluctuations ^Pa(b) should each contain a half of the origi- 
nal vortices in and a half of antivortices. This is readily 
achieved by including those (anti)vortex variables whose 
positions are labeled by with i even into {pA, while the 
odd ones are absorbed into lpb- The symmetrization is 
just a convenient mathematical tool that automatically 



guarantees this goal. 

The above symmetrization leads to the new partition 
function: Z — s- Z: 

Z ^ j V¥ j j 25^^ j exp [- dr J d^rC] , 

(9) 

in which the half- flux-to-minus- half- flux (Z2) symmetry 
of the singular gauge transformation is now manifest: 

C = ¥[{dr + iar)cro + ivras]^ + ^^W* + Co[Vf„ a^], 

(10) 

where Cq is the "Jacobian" of the transformation given 
by 

e-fo'drfd-rco ^2-'''J2 I T^^ir.r) (11) 



Here are the Pauli matrices, (3 = 1/T, H' is given in 
Eq. (P) and, for later convenience, £0 will also include 
the energetics of vortex core overlap driven by amplitude 
fluctuations and thus independent of long range super- 
flow (and of A). We call the transformed quasiparticles 
»I/t — (-01,^1) appearing in ( [l0| ) "topological fermions" 
(TF). TF are the natural fermionic excitations of the 
pseudogapped normal state. They are electrically neu- 
tral and are related to the original quasiparticles by the 
inversion of transformation (^. 

By recasting the original problem in terms of topologi- 
cal fermions we have accomplished our original goal: the 
interactions between quasiparticles and vortices are now 
described solely in terms of two local superflow fields [E3| : 



VAi^L = d^ipA 



(12) 



which we can think of as superfluid velocities associated 
with the phase configurations (Pa{b){^,t) of a (2-1-1)- 
dimensional XY model with periodic boundary condi- 
tions along the r-axis. Our Doppler and Berry gauge 
fields and are linear combinations of VAf^ and . 
Note that is produced exclusively by vortex defects 
since the "spin-wave" configurations of ip can be fully 
absorbed into . All that remains is to perform the sum 
in (||) over all the "spin-wave" fluctuations and all the 
spacetime configurations of vortex defects {ri(r)} of this 
(2-|-l)-dimensional XY model. 



C. "Coarse-grained" Doppler and Berry U(l) 
gauge fields (ii^ and a^) and their physical 
significance 

Unfortunately, the exact integration over the phase 
Lp{r,T) is prohibitively difficult. To proceed by analytic 
means we must devise some approximate procedure to 
integrate over vortex-antivortex positions {ri{T)} in (|^) 
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which will capture qualitative features of at least the 
long-distance, low-energy physics of the original prob- 
lem. This is where our recasting of the problem in terms 
of BdG fermions interacting with superflow fields VAfi 
and VBfj. ( |T^ ) will come in handy. A hint of how to de- 
vise such an approximation comes from examining the 
role of the Doppler gauge field in the physics of this 
problem. To illustrate our reasoning and for simplicity, 
we consider the finite T case where we can ignore the 
T-dependence of (p(r,T). The results are easily general- 
ized, with appropriate modifications, to include quantum 
fluctuations. 

We start by noting that Vj = 2v — (2e/c)A is just 
the physical superfluid velocity ^ , invariant under both 
A ^ B singular gauge transformations (^) and ordi- 
nary electromagnetic U(l) gauge symmetry. The su- 
perfluid velocity (Doppler) field, swirling around each 
(anti)vortex defect, is responsible for a vast majority of 
phenomena that we associate with vortices: long range 
interactions between vortex defects, coupling to exter- 
nal magnetic field and the Abrikosov lattice of the mixed 
state, Kosterlitz-Thouless transition, etc. Remarkably, 
we will show that its role is essential even for the physics 
discussed in the present paper, although it now appears 
as a supporting actor to the Berry gauge field which 
ultimately occupies the center stage. To see how this 
comes about, imagine that a (a^) were absent - then, 
upon integration over topological fermions in (|l0|), we 
obtain the following term in the effective Largangian: 

M^I^V^-^A^ + (...) , (13) 

where (■ • • ) denotes higher order powers and deriva- 
tives of — (2e/c)A. In the above we have replaced 
V (1/2) V93 to emphasize that the leading term, with 
the coefficient A'P proportional to the bare superfluid 
density, is just the standard superfluid-velocity-squared 
term of the continuum XY model - the notation for 
the coefficient will become clear in a moment. We can 
now write Vi^ — Viy9vortcx + Viy9spin-wavc and reexpress 
the transverse portion of (|l3|) in terms of (anti)vortex 
positions {r^} to obtain a familiar form: 




or, by using V • V(y9vortex = and V x Vi^vortox = 27rp(r), 
equivalently as: 

J d\ J d^r'p{r)p{r')\n\r-r'\ , (15) 

where p{r) = 27r^j qiS{r — r^) is the vortex density. 

The Meissner coupling of v to fermions is very strong - 
it leads to familiar long range interactions between vor- 
tices which constrain vortex fluctuations to a remark- 
able degree. To make this statement mathematically ex- 
plicit we introduce the Fourier transform of vortex den- 
sity p(q) = J2i exp(zq • r^) and observe that its variance 



satisfles: 

(p(q)p(-q)) « ■ (16) 

Vortex defects form an incompressible liquid - the long 
distance vorticity fluctuations are strongly suppressed. 
This "incompressibility constraint" is naturally enforced 
in (|^) by replacing the integral over discrete vortex posi- 
tions {rj} with the integral over continuously distributed 
field p(r) with (p(r)) — 0. The Kosterlitz-Thouless tran- 
sition and other vortex phenomenology are still main- 
tained in the non-trivial structure of £o[p(r)]. But the 
long wavelength form of (^3|) now reads: 

M2 (^2v-^A^ +(•••) , (17) 

and can be interpreted as a massive action for a U(l) 
gauge field v. The latter is our coarse grained Doppler 
gauge field defined by 

V X v(r) = 7rp(r) . (18) 

We have now gone full circle with the Doppler gauge field 
V. The coarse-graining procedure has made it into a mas- 
sive U(l) gauge field v whose influence on TF disappears 
in the long wavelength limit. We can therefore drop it 
from low energy fermiology. Hereafter we shall drop the 
bar and use symbol v for both the actual and the coarse 
grained quantity, the meaning being obvious from the 
context. 

The real problem also contains the Berry gauge fleld 
a (op) which now must be restored. However, if we are to 
take advantage of introducing continuous p(r) and even- 
tually dispensing with v, a (a^) cannot remain the same 
gauge field we started with in Eq. (^^ when we embarked 
on our quest to derive the effective theory. Having re- 
placed the Doppler field v with the distributed quantity 
(p8|), we cannot simply continue to keep a (a^J speci- 
fied by half-fiux Dirac (Aharonov-Bohm) strings located 
at discrete vortex positions {r^}. Instead, a (a^) must 
be replaced by some new gauge field which refiects the 
"coarse-graining" that has been applied to v - simply 
put, the Z2- valued Berry gauge field ^ of the original 
problem (^0|) must be "dressed" in such a way as to best 
compensate for the error introduced by "coarse-graining" 
V. In the language of RG, we must find such "dressing" 
of the Berry gauge field, i.e. the form and the bare action 
for a (flp), which renders any such errors unimportant for 
the low energy physics. 

The recipe for such a required "dressing" of the Berry 
gauge field is straightforward in the FT singular gauge 
- it takes the form of a non-compact U(l) gauge field 
with a simple Maxwellian action. To see how this comes 
about note that if we insist on replacing v by its "coarse- 
grained" form (|l^), the only way to achieve this is to 
"coarse-grain" both va and in the same manner: 

V X ^ 2npA ; V X VB ^ 2npB , (19) 
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where ^^^^^(r) are now continuously distributed densi- 
ties of A{B) vortex defects (the reader should contrast 
this with This is because the elementary vortex 

variables of our problem are not the sources of v and a; 
rather, they are the sources of and v^. We cannot 
separately fluctuate or "coarse-grain" the Doppler and 
the Berry "halves" of a given vortex defect - they are 
permanently confined into a physical {hc/2e) vortex. On 
the other hand, we can independently fluctuate A and B 
vortices - this is why it was important to use the singular 
gauge (^) to rewrite the original problem solely in terms 
of A and B vortices and associated superflow fields (12). 

Following this recipe we can now reassemble the coarse- 
grained Doppler and Berry gauge fields as: 



1/ ^ 1/ ^ 

with the straightforward generalization to (2+l)D: 



:{vA - vb) 



(20) 



(21) 



The coupling of "vaIva) and vb(wb) to fermions is a 
hybrid of Meissner and minimal coupling [p7| . They 
contribute a product of U(l) phase factors exp(i(5^) x 
exp{i6B) to the BdG fermion loops with both exp(iSA) 
and exp(iSB) being "random" in the sense of previous 
subsection. Upon coarse-graining va(va) and vb{vb) 
turn into non-compact U(l) gauge fields and therefore 
v(t;) and a(a) must as well. 



D. Further remarks on FT gauge 

The above "coarse-grained" theory must have the fol- 
lowing symmetry: it has to be invariant under the ex- 
change of spin- up and spin-down labels, ipi^ ^ ip^, with- 
out any changes in v (v). This symmetry insures that 
(arbitrarily preselected) Sz component of the spin, which 
is the same for TF and real electrons, decouples from the 
physical superfluid velocity which naturally must couple 
only to charge. When dealing with discrete vortex de- 
fects this symmetry is guaranteed by the singular gauge 
symmetry defined by the family of transformations (|^). 
However, if we replace and vs by their distributed 
"coarse-grained" versions, the said symmetry is preserved 
only in the FT gauge. This is seen by considering the ef- 
fective Lagrangian expressed in terms of coarse-grained 
quantities: 



C ^ ¥dr^ + ¥H'i> + Co[pA,PB] 



(22) 



where 7i' is given by (^, v^, vs are connected to pA, 
Pb via (p^, and Cq is independent of r. 

The problem lurks in £o[pA,Ps] ^ this is just the 
entropy functional of fluctuating free A{B) vortex- 
antivortex defects and has the following symmetry: 
Pa{b) ~^ ~Pa(b) with Pb{a) kept unchanged. This sym- 
metry reflects the fact that the entropic "interactions" do 



not depend on vorticity. Above the Kosterlitz-Thouless 
transition we can expand: 

£0 - ^(V X v^)2 + :^(V X vs)^ + (•••), (23) 



where the ellipsis denote higher order terms and the co- 
efficients K^^g^ "f*'^'' (see Appendix A for details). 

The above discussed symmetry of Hamiltonian Ti.' (|^) 
demands that v and a (|2^,|2^) be the natural choice for 
independent distributed vortex fluctuation gauge fields 
which should appear in our ultimate effective theory. Cq, 
however, collides with this symmetry of Ti.' - if we replace 
'^A(B) ^ v±a in (^3|), we realize that v and a are coupled 
through Cq in the general case Ka 7^ Kb- 



Co 



Ka + Kb 



((V X v)2 + (V X a)2) 



+ ^£_^(V X a) . (V X V) 



(24) 



Therefore, via its coupling to a, the superfluid velocity 

V couples to the "spin" of topological fermions and ul- 
timately to the true spin of the real electrons. This is 
an unacceptable feature for the effective theory and se- 
riously handicaps the general "A — i?" gauge, in which 
the original phase is split into cp ipA + fB, with (pA(B) 
each containing some arbitrary fraction of the original 
vortex defects. In contrast, the symmetrized transforma- 
tion (^ which apportions vortex defects equally between 
ipA and (fiB ||2^ leads to Ka = Kb and to decoupling of 

V and a at quadratic order, thus eliminating the prob- 
lem at its root. Furthermore, even if we start with the 
general "A — i?" gauge, the renormalization of Cq aris- 
ing from integration over fermions will ultimately drive 
Ka — Kb — *■ and make the coupling of v and a irrele- 
vant in the RG sense |2^ . This argument is actually quite 
rigorous in the case of quantum fluctuations where the 
symmetrized gauge (E ) represents a fixed point in the RG 
analysis (see below) 1 30| . Consequently, it appears that 
the symmetrized singular transformation (^) employed in 
( p^ ) is the preferred gauge for the construction of the ef- 
fective low-energy theory |3^. In this respect, while 
all the singular A — B gauges are created equal some are 
ultimately more equal than others. 

The above discussion provides the rationale behind us- 
ing the FT transformation in our quest for the effective 
theory. At low energies, the interactions between quasi- 
particles and vortices are represented by two U(l) gauge 
fields V and a (|20| , ^ ). The conversion of a from a Z2- 
valued to a non-compact U(l) field with Maxwellian ac- 
tion is effected by the confinement of the Doppler to the 
Berry half of a singly quantized vortex ~ in the coarse- 
graining process the phase factors exjp{iS) of the non- 
compact Doppler part "contaminate" the original (±1) 
factors supplied by a (|^) . This contamination diminishes 
as doping a; — > since then vp/vA and "vortices" 
are effectively liberated of their Doppler content. In this 
limit, the pure Z2 nature of a is recovered and one enters 



9 



the realm of the 1i gauge theory of Senthil and Fisher 
p2| . In contrast, in the pairing pseudogap regime of this 
paper where wf/^a > 1 and singly quantized (hc/2e) 
vortices appear to be the relevant excitations |33| , we ex- 
pect the effective theory to take the U(l) form described 
by V and a 



E. Jacobian Co[Vfj,,afi] and its Maxwellian form 

Having elucidated the origin of the "coarse-grained" 
U(l) gauge fields and and settled on the sym- 
metrized FT transformation as the natural gauge 
choice for this problem, there remains one more task 
to be accomplished before we can conclude this Section. 
We need to derive a precise expression for the long dis- 
tance, low energy form of the "Jacobian" Co[Vfj,, a^] which 
serves as the "bare action" for the gauge fields and 
of our effective theory. As shown below, this form is a 
non-compact Maxwellian whose stiffness K (or inverse 
"charge" — K) stands in intimate relation to the 

helicity modulus tensor of a dSC and, in the pseudogap 
regime of strong superconducting fluctuations, can be ex- 
pressed in terms of a finite physical superconducting cor- 
relation length Csc, K (X (2D); K oc ((2-hl)D). 
As we enter the superconducting phase and 4c 00, 
K 00 ds. well (or — *■ 0), implying that and 
have become massive. Our derivation, the results 
of which were originally quoted and used in Ref. [p^ , 
can be accomplished with remarkably little algebra and 
holds for Ginzburg-Landau, XY model or any other rep- 
resentation of superconducting fluctuations. This is no 
accident - the straightforward relationship between the 
massless (or massive) character of Lo[v^,a^] and the su- 
perconducting phase disorder (or order) is a consequence 
of rather general physical and symmetry principles. 

To make good on the above claim consider first a sim- 
ple example of an s-wave superconductor with a large 
gap A extending over all of the Fermi surface. We can 
also view this as a model for the high energy BdG quasi- 
particles in dSC, those far removed from the nodes. The 
action takes the form similar to (p^: 



but with 7i(, defined as: 



(25) 



(26) 



with TT — f>-\-aL. in the fermionic action (but not in £0) 
goes into Vf^ — {e/c)Af^ when the electromagnetic field is 
included. The reader might wish to recall here that we 
have defined Co[vij,,aij] in this particular way to clearly 
separate the superflow mediated interactions among vor- 
tices, which include A^, from the entropic effects and 
short range amplitude driven core-overlap interactions, 
which do not. 



1. Thermal vortex-antivortex fluctuations 

We can now reap the benefits of this convenient sep- 
aration. In the language of BdG fermions the system 
( p5| ) is a large gap "semiconductor" and the Berry gauge 
field couples to it minimally through "BdG" vector and 
scalar potentials a and a^-. Such BdG semiconductor is 
a poor dielectric diamagnet with respect to a^. We pro- 
ceed to ignore its "diamagnetic susceptibility" and also 
set Or = to concentrate on thermal fluctuations. All 
this means is that the Berry gauge field part of the cou- 
pling between quasiparticles and vortices in a large gap 
s-wave superconductor influences the latter only through 
weak short-range interactions which are unimportant in 
the region of strong vortex fluctuations near Kosterlitz- 
Thouless transition. We can therefore drop a from the 
fermionic part of the action and integrate over it to ob- 
tain £0 in terms of physical vorticity p{r) = (V x v)/7r. 
Additional integration over the fermions produces the ef- 
fective free energy functional for vortices: 



c 



r + {■■■) +^o[p] 



(27) 



where we have used our earlier notation and have intro- 
duced a small external transverse vector potential A'^'^*. 
The ellipsis denotes higher order contributions to vor- 
tex interactions. As discussed earlier in this Section, 
the familiar long range interactions between vortices lead 
directly to the standard Coulomb gas representation of 
the vortex-antivortex fluctuation problem and Kosterlitz- 
Thouless transition. 

The presence of these long range interactions im- 
plies that the vortex system is incompressible ( p^ ) and 
long distance vortex density fluctuations are suppressed. 
When studying the coupling of BdG quasiparticles to 
these fluctuations it therefore suffices to expand the "en- 
tropic" part: 



X v)2 + 



(28) 



The above expansion is justified above Tc since we know 
that at T 3> Tc we must match the purely entropic form 
of a non-interacting particle system, Cq oc Sp^ (see Ap- 
pendix A). 

To uncover the physical meaning of coefficient K we 
expand the free energy F of the vortex system to second 
order in A'^'^*. In the pseudogap state the gauge invari- 
ance demands that F depends only on V x A"'^': 



F[V X A° 



(2e)! 
2c2 



X I d^r{V X A°''*)2 



(29) 

Note that {{2e)'^ / c'^)x is just the diamagnetic susceptibil- 
ity in the pseudogap state, x determines the long wave- 
length form of the helicity modulus tensor T^i,(g) defined 



S^F 



,5A-t(9)<5^"/H-9) 



(30) 
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The above is the more general form of T^,y (q) appUcable 
to uniaxially symmetric 3D and (2+l)D XY or Ginzburg- 
Landau models; in 2D only T^y(q) appears. In the long 
wavelength limit T ^,^{q) vanishes as q^X'- 



i2ey 



(31) 



for the isotropic case, while for the anisotropic situation 

X± 7^ X||: 

^2 

rT^i/(q) = (x|l - X±)^zaf,ezi3uqaqp + X±^pap.<^P0^qaqfj 

(32) 



eajS'Y is the Levi-Civita symbol, summation over repeated 
indices is understood and index z of the anisotropic 3D 
GL or XY model is replaced by r for the (2+l)D case. 
(2e)^/c^ is factored out for later convenience. 

So, what is K7 Let us compute the helicity modulus 
of the problem explicitly. This is done by absorbing the 
small transverse vector potential A'''^* into v in (|2^,^) and 
integrating over the new variable v — (e/c)A°'^*. The he- 
cility modulus tensor measures the screening properties 
of the vortex system. In a superconductor, with topologi- 
cal defects bound in vortex-antivortex dipoles, there is no 
screening at long distances. This translates into Meissner 
effect for A®'^*. When the dipoles unbind and some free 
vortex-antivortex excitations appear the screening is now 
possible over all lengthscales and there is no Meissner ef- 
fect for A®''*. The information on presence or absence 
of such screening is actually stored entirely in Cq^ where 
^ext reemerges after the above change of variables. We 
finally obtain: 



Ax = K 



T 



lim 

|qHO+ 



dVe'i'"((5/9(r)(5p(0)) , (33) 



where the thermal average (• • • ) is over the free energy 
( p7| ) with A°''* = or, equivalently, (|25|). In the normal 
phase only the first term contributes in the long wave- 
length limit, the second being down by an extra power 
of q^ courtesy of long range vortex interactions. Conse- 
quently, K — Ax (the factor of 4 is due to the fact that 
the true superfluid velocity is 2v [p8|). 

We see that in the pseudogap phase, with free vor- 
tex defects available to screen, Co[v] takes on a massless 
Maxwellian form, the stiffness of which is given by the 
diamagnetic susceptibility of a strongly fluctuating su- 
perconductor. In the fluctuation region x is given by the 
superconducting correlation length p4|: 



(34) 



where C2 is a numerical constant, intrinsic to a 2D GL, 
XY or some other model of superconducting fluctuations. 

As we approach Tc, ~f 00 and the stiffness of 
Maxwellian term ( p^ ) diverges. This can be interpreted 
as the Doppler gauge field becoming massive. Indeed, 



immediately bellow the Kosterlitz-Thouless transition at 
Tc, £0 — ™^v^ -I- . . . , where m„ ^ M and m„ — > as 
T . This is just a reflection of the helicity modulus 

tensor now becoming finite in the long wavelength limit, 

46^ ■mlM'^ 



c2 -f ml ' 



Topological defects are now bound in vortex-antivortex 
pairs and cannot screen resulting in the Meissner effect 
for A'"'*. The system is a superconductor and v had 
become massive. 

Returning to a d-wave superconductor, we can re- 
trace the steps in the above analysis but we must re- 
place TL'g in ( p5| ) with TL' (||). Now, instead of a large 
gap BdG "semiconductor" , we are dealing with a nar- 
row gap "semiconductor" or BdG "semimetal" because 
of the low energy nodal quasiparticles. This means that 
we must restore the Berry gauge field a to the fermionic 
action since the contribution from nodal quasiparticles 
makes its BdG "diamagnetic susceptibility" very large, 
XBdG ~ l/T > l/T* (see the next Section). The long 
distance fluctuations of v and a are now both strongly 
suppressed, the former through incompressibility of the 
vortex system and the latter through XBdG- This allows 
us to expand (k3): 



(35) 



/:o^^(Vxv^)^ + :^(Vxv5)2 + (...) 



where Ka = Kb = K \s mandated by the FT singular 
gauge. Since in our gauge the fermion spin and charge 
channels decouple A'^'^* still couples only to v and the 
above arguments connecting K to the helicity modulus 
and diamagnetic susceptibility x follow through. This 
flnally gives the Maxwellian form of Rcf . |13] : 



A-f (Vxv)2 + |(Vxa)^ 



(36) 



where K is still given by Eq. (pi]). Note however that 
^sc of a d-wave and of an s-wave superconductor are two 
rather different functions of T, x and other parameters of 
the problem, due to strong Berry gauge field renormaliza- 
tions of vortex interactions in the d-wave case. Nonethe- 
less, as T ^ Tc, the Kosterlitz-Thouless critical behavior 
remains unaffected since XBdGi while large, is still finite 
at all finite T . 

Just as advertised, we have shown that Cq[pat Pb\ — 
£o[v, a] in the pseudogap state takes on the massless 
Maxwellian form (^^, with the stiffness K set by the 
true superconducting correlation length £_sc (|3^). This 
result holds as a general feature of our theory irrespec- 
tive of whether one employs a Ginzburg-Landau the- 
ory, XY model, vortex-antivortex Coulomb plasma, or 
any other description of strongly fluctuating dSC, as 
long as such description properly takes into account 
vortex-antivortex fluctuations and reproduces Kosterlitz- 
Thouless phenomenology. In the Appendix A we show 
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that within the continuum vortex-antivortex Coulomb 
plasma model: 



T 



27r2 



[(V X v)2 + (V X a)2] , 



(37) 



where n; is the average density oi free vortex and antivor- 
tex defects. Comparison with (36) allows us to identify 



2. Quantum fluctuations of (2+l)D vortex loops 

The above results can be generalized to quantum fluc- 
tuations of spacetime vortex loops. The superflow fields 



VA(B)t, satisfy {d x va^i 



where jA(i3)p(a;) 



are the coarse-grained vorticities associated with A{B) 
vortex defects and = 0. The topology of vor- 

tex loops dictates that jA{B)tiix) be a purely transverse 
field, i.e. d ■ jA(B) — 0; reflecting the fact that loops 
have no starting nor ending point. Again, we begin 
with a large gap s-wave superconductor and use its poor 
BdG diamagnetic/dielectric nature to justify drop ping 
the Berry gauge field from the fermionic part of (|25| ) 
and integrating over in the "entropic" part containing 

The integration over the fcrmions contains an impor- 
tant novelty specific to the (2-|-l)D case: the appearance 
of the Berry phase terms for quantum vortices as they 
wind around fermions. Such Berry phase is the conse- 
quence of the first order time derivative in the original 
fermionic action ([^). If we think of spacetime vortex 
loops as worldlines of some relativistic quantum bosons 
dual to the Cooper pair field A(r, t), as we do in the Ap- 
pendix A, then these bosons see Cooper pairs and quasi- 
particles as sources of "magnetic" flux . At the mean 
field level, this translates into a dual "Abrikosov lattice" 
or a Wigner crystal of holes in a dual superfluid. Ac- 
cordingly, the non-superconducting ground state in the 
pseudogap regime will likely contain a weak charge modu- 
lation - the modulation is made weak by the same strong 
fluctuations that make Tc <C r*(rNernst)- The focus of 
the present paper being the symmetric AFL description 
of the pseudogap, we postpone the discussion of this point 
to Part II and will ignore it for the rest of this paper. 
This is justified by the fact that does not couple to 
charge directly and is quantitatively valid in the window 
T, A/ ^ (a;,WF,A|q|) <C T* , where A/ is any small gap 
in the nodal TF spectrum produced by the said charge 
modulation. 

Hereafter, we blissfully turn the blind eye to the above 
subtleties and assume that the transition from a dSC into 
a pseudogap phase proceeds via unbinding of vortex loops 
of a (2-|-l)D XY model or its GL counterpart or, equiva- 
lently, an anisotropic 3D XY or GL model where the role 
of imaginary time is taken on by a third spatial axis z 
psf . Having learned all we really need from an earlier 2D 
example we can now integrate the fermions to obtain the 



effective Lagrangian for coarse-grained spacetime loops 
(ttj^ = (a x v)^): 

C[j,] = ACi2v, -Al^r + (•••) + Colj,] , (38) 



where = My = M and = M/cs, with Cg ^ 
vp being the effective "speed of light" in the vortex 
loop spacetime. The incompressibility condition reads 
{jt^{q)ji'i-<l)) ~ q'^i^f^i' - {qtJ.qv I ff')] and in the pseudo- 
gap state permits the expansion: 



(39) 



where — Ky = K ^ . Using the analogy with the 
uniaxially symmetric anisotropic 3D XY (or GL) model 
we can expand the ground state energy in the manner of 



(2e)^ 
2c2 

_L, 



E{d X = £;[0] 



X_L.li / <rx{d X A^'^y^. 



(40) 



with x_L = Xx =_Xy = X and X|| ^ Xt ^ X- Note that 
the form of Lq ( p9|) follows directly from the requirement 
that there are infinitely large vortex loops, resulting in 
vorticity fluctuations over all distances. Combined with 
( ^0| ) and then translated to the language of a (2-|-l)D XY 
(GL) model it tells us something already familiar: upon 
the transition to the pseudogap state generated by vor- 
tex loops unbinding, the superconductor has turned into 
an insulator [Q. x and Xt determine the diamagnetic 
and dielectric susceptibilities of this insulating pseudogap 
state. 

The explicit computation of T^^{q) ( |3C| , |3^ ) leads to: 



4x^ 



K,r 



lim 



(41) 



where the second term is again eliminated by the incom- 
pressibility of the vortex system. This results in: 



if = 4x = 4C3er ; Xr = 4xr = 4C: 



(42) 



where we used the result for the anisotropic 3D XY 
or GL models: Xi = CaTe,,, X|| = C^Xejiw. Cg be- 
ing the intrinsic numerical constant for those models 
(C3 7^ C2). In the case of (2+l)D vortex loops ^t- cx ^sc 
since our adopted model has the dynamical critical ex- 
ponent z = 1. Thus, we again encounter a massless 
Maxwellian form ( ^9[ ) whose stiffness diverges as we ap- 
proach a superconductor except now this divergence is 
linear in the superconducting correlation length, K S^^c- 
The application to a d-wave superconductor is straight- 
forward: the nodal structure of BdG quasiparticles in 
dSC helps along by the way of providing the anomalous 
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stiffness for the Berry gauge field ~ upon integration 
over the nodal fcrmions the following term emerges in the 
effective action: 

^ \<l\[5^L^ - {q^iq,^lq^)]ai,{q)a^{-q) ■ (43) 

In the terminology of our chimerical BdG "semimetal" , 
the "susceptibility" XBdG is not merely very large, it di- 
verges: XBdG ^ l/^i as g — > 0, and is computed in detail 
in Section IV. We use this pleasing fact to observe that 
we are fully justified in expanding ColjA^niB^i] and re- 
taining only quadratic terms as long as we keep a safe 
distance from the actual phase transition: 

C,^^{d^VA)l + ^{d^vs)l + {---) , (44) 

where Ka^, = Kb^ = K^^ is again assured by our choice 
of the FT singular gauge ^^). 

The above reasoning merits an amusing aside: beside 
the ubiquitous incompressibility {{{d x v)fj,{d x v)„) ~ 
q'^i^fj.f ~ {qfj.qiy/q'^)]) the integration over nodal fermions 
now also occasions diverging xbcIG implying ((9x a)^(9x 
a)i/) ~ q[^^ll' — [qtiqv/q^)]- Had we chosen a singular 
gauge in which the sets A and B were not equivalent, 
like Anderson gauge and therefore K^n^ Kb^, the 
ensuing {dxv)-{dx a) coupling in Co ( p4| , p4[ ) would now 
be driven to zero in the long distance limit as an irrelevant 
operator in the RG sense ||3q| . The reason is simple: the 
coupling of gauge fields v^j. and to topological fermions 
mandates that they decouple at the quadratic (harmonic) 
level due to decoupling of "charge" and "spin" channels 
for TF. Such coupling of and can only arise from 
Cq by our uninformed choice of a singular gauge. Since 
both (9 X u)^ and {d x a)^ terms in ( ^^ are strongly 
relevant due to diverging contributions they receive from 
fermions, the coupling constant in front of {d x v) ■ {d x 
a), proportional to Kai^ ~ Kb^j., is driven to zero under 
repeated applications of RG transformation. Therefore, 
the FT gauge (^|l]) specifically designed to insure Kaii — 
Kbii = at the very start, is recovered as an RG fixed 
point 

Finally, we rewrite (^) in terms of , ~ 
(l/2)(t;AM±t>s^): 

C,^^{dxv)l + ^{dxa)l , (45) 

and observe that the fact that A'j^^ couples only to 
means that the expression (^) for is still valid. Of 
course, ^sc(2;,T) is now truly different from its s-wave 
counterpart, including a possible difference in the critical 
exponent, since the coupling of d-wave quasiparticles to 
the Berry gauge field is marginal at the RG engineering 
level and may change the quantum critical behavior of 
the superconductor-pseudogap (insulator) transition. 



FIG. 2: Schematic representation of the Fermi surface of the 
cuprate superconductors with the indicated nodal points of 
the d^2_y2 gap. 



III. QEDa ~ A LOW ENERGY EFFECTIVE 
THEORY OF THE PSEUDOGAP STATE 

We have now elucidated the nature of the coupling of 
our two gauge fields, and a^, to TF and have specified 
their "bare" thermal and quantum dynamics encoded in 
Co [Vfj^ , a^] (|3^,^) . To make further progress toward our 
ultimate goal of describing the low energy fermiology in 
the pseudogap state, we now focus our attention on the 
nodal quasiparticle excitations of the Hamiltonian Ti.' Eq. 

. This will enable us to apply the machinery of pertur- 
bative RG to nodal (massless) TF and rid our effective 
theory of all remaining excess baggage. 



A. Farewell to and residual interactions 

As indicated in Fig. H, the low energy quasiparticles are 
located at the four nodal points of the dxy gap function: 

I = {±kF,0) and k2.2 — (0,±/ci?), hereafter denoted 
as (1,1) and (2,2) respectively. To focus on the leading 
low energy behavior of the fermionic excitations near the 
nodes we follow the standard procedure |]37|| and linearize 
the Lagrangian (|l^). To this end we write our TF spinor 
^ as a sum of 4 nodal fermi fields, 

(46) 

The {72 matrices have been inserted here for convenience: 
they insure that we eventually recover the conventional 
form of the QED3 Lagrangian. (Without the (T2 matri- 
ces the Dirac velocities at 1,2 nodes would have been 
negative.) Inserting 'F into ( p^ ) and systematically ne- 
glecting the higher order derivatives p^, we obtain the 
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nodal Lagrangian of the form 



C-D = ^¥^{Dr^ivFD^aj,+iv,^Dya^\^i (47) 



1=2,2 



where = 9^ 
and £o [a/j] isgiven by (fl5| 



ia^ denotes the covariant derivative 



ii^^(axa)2^-i^(5xa)2 



(48) 



vf = dek/dk denotes the Fermi velocity at the node and 
va = dAi^/dk denotes the gap velocity. Note that vp 
and va already contain renormalizations coming from 
high energy interactions and are effective material pa- 
rameters of our theory. Similarly, if^ = 1/e^, derived in 
the previous section in terms of ^sc(a:^, T), are treated as 
adjustable parameters which are matched to experimen- 
tally available information on the range of superconduct- 
ing correlations in the pseudogap state. 

Doppler gauge field has disappeared from the above 
expression. After informing us on how to properly 
"coarse-grain" the theory and dressing our Z2-valued 
Berry gauge field in its ultimate U(l) Maxwellian outfit, 
the time has come to drop w^, its eventual demise caused 
by the Meissner coupling to BdG fermions discussed in 
the previous section. After being "screened" by high en- 
ergy and nodal TFs it is rendered massive both in the 
superconducting and pseudogap states and unimportant 
for low energy physics. Its legacy lives on, however, hav- 
ing given birth to the U(l) non-compact character of a^. 

In contrast. Berry gauge field remains massless in 
the pseudogap state, as it cannot acquire mass by cou- 



pling to the fermions. As seen from Eq. 



Cu couples 



minimally to the Dirac fermions and therefore its mass- 
less character is protected by gauge invariance. Physi- 
cally, one can also argue that couples to the TF spin 
three-current - in a spin-singlet d-wave superconductor 
SU(2) spin symmetry must remain unbroken, thereby in- 
suring that On remains massless. Its massless Maxwellian 
dynamics (Eq) in the pseudogap state can therefore be 
traced back to the topological state of spacetime vortex 
loops and directly reflect the absence of true supercon- 
ducting order (Section II) or, equivalently, the presence 
of "vortex loop condensate" and dual order (Appendix 
A). 

We have also dispensed with the residual interactions 
represented by Tires in (§) . These interactions are generi- 
cally short-ranged contributions from the p-h and ampli- 
tude fluctuations part of the p-p channel and in our new 
notation are exemplified by: 



Hrc 



(49) 



1,1' 



The effective vertex Iw has a scaling dimension —1 at 
the engineering level. This follows from the RG analysis 



near the massless Dirac points which sets the dimension 
of to [length]"^. The implication is that is ir- 
relevant for low energy physics in the perturbative RG 
sense and we are therefore well within our rights in set- 
ting ^ 0. However, the residual interactions will not 
go so quietly into the night: such interactions are known 
to become important if stronger than some critical value 
Ic . In the present theory, this is bound to happen 
in the severely underdoped regime and at half-filling, as 
X ^ ll^. In this case, the residual interactions are 
becoming large and comparable in scale to the pairing 
pseudogap A, and are likely to cause chiral symmetry 
breaking (CSB) which leads to spontaneous mass gen- 
eration for massless Dirac fermions. The CSB and its 
variety of patterns in the context of the theory (47 4^) 
in underdoped cuprates were discussed in ||^, ^2[ and 
are the subject of Part II of this paper. Here, where we 
have limited ourselves to the chirally symmetric phase, 
we assume < Ic, which we expect to be the case for 
moderate underdoping, and set 0. 

In the end, we are left with (|^,||) as our effective low 
energy theory for nodal TF. This theory, derived pre- 
viously in Ref. fl^ , is the chief dynamical muscle be- 
hind the physics discussed in this paper. It describes the 
problem of massless topological fermions interacting with 
massless vortex "Beryons" , i.e. the quanta of the Berry 
gauge field a^, and is formally equivalent to the Euclidean 
quantum electrodynamics of massless Dirac fermions in 
(2-fl) dimensions (QED3). It, however, suffers from an 
intrinsic Dirac anisotropy by virtue of Vi;- ^ va- 



B. QED3 Lagrangian for the pseudogap state 

We are now in position to do some real calculations 
within our theory. Before we plunge into the algebra, 
however, we first apply some cosmetics: Lagrangian ( |47| ) 
is not in the standard from as used in quantum electro- 
dynamics where the matrices associated with the com- 
ponents of covariant derivatives form a Dirac algebra 
and mutually anticommutc. In ( p7| ) the temporal deriva- 
tive is associated with unit matrix and it therefore com- 
mutes (rather than anticommutes) with ai and (73 ma- 
trices associated with the spatial derivatives. These non- 
standard commutation relations, however, lead to some 
rather unwieldy algebra. For this reason we manipulate 
the Lagrangian ( p7| ) into a slightly different form that 
is consistent with usual field-theoretic notation. First, 
we combine each pair of antipodal (time reversed) two- 
component spinors into one 4-component spinor. 



T, = 



*2 
*2 



(50) 



Second, we define a new adjoint four-component spinor 
= -zT|7o. (51) 
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In terms of this new spinor the Lagrangian becomes 

i^-D^Y. ^n.Dl^^r, + ^K^{d X a)l , (52) 

1=1,2 

with covariant derivatives 

D^^^ = i[{dr + iar),VFidx +iax),VA{dy +iay)], 
-D^f' = i[idr + iar),VFidy + iay),VAidx + iax)]- 
The 4x4 gamma matrices, defined as 



70 = 

72 = 



0-2 

-^2 

-as 
as 



71 = 



CTl 







-CTl 



now form the usual Dirac algebra, 

{7m,7i/} = 2(5^1^ 
and furthermore satisfy 

Tr(7^) = , Tr(7^7,) = 4<5^,. 



(53) 



(54) 



(55) 



The use of the adjoint spinor T instead of conven- 
tional is a purely formal device which will simplify 
calculations but does not alter the physical content of 
the theory. At the end of the calculation we have to re- 
member to undo the transformation ( [5l| ) by multiplying 
the (T(x)T(x')) correlator by ijo to obtain the physical 
correlator {T{x)T\x')). 

Next, to make the formalism simpler still we can elim- 
inate the asymmetry between the two pairs of nodes by 
performing an internal SU(2) rotation at nodes 2, 2: 



T2 ^ e-^-^°7iT2, 
leading to the anisotropic QED3 Lagrangian 



Cd=J2 - + ^K^id X a)l 

1=1,2 



(56) 



(57) 



IV. SPECTRAL PROPERTIES OF 
TOPOLOGICAL FERMIONS AND PHYSICAL 
ELECTRONS IN QED3 

We shall start by considering the isotropic case, vp = 
t!A = 1, which although unphysical in the strictest sense, 
is computationally much simpler and provides penetrat- 
ing insights into the physics embodied by the QED3 La- 
grangian (^). After we have understood the isotropic 
case we will then be ready to tackle the calculation for 
the general case and will show that Dirac cone anisotropy 
does not modify the essential physics discussed here. To 



a) 



q + k 




b) 



k — q 



FIG 
(b). 



3: One loop Berryon polarization (a) and TF self energy 



make contact with standard literature on QED3, we fur- 
ther consider a more general problem with N pairs of 
nodes described by Lagrangian 



N 

E 



^Ht^iidf, ~ a^)T/ + ^Kf,{d x a) 



(58) 



For the basic problem of a single Cu02 layer N = 2. 
As we will show in the next section, N itself is variable 
and can be equal to four or six in bi- and multi-layer 
cuprates. Our analytic results can be viewed as arising 
from the formal expansion, although we expect them 
to be qualitatively (and even quantitatively!) accurate 
even for A'^ = 2 as long as we are within the symmetric 
phase of QED3 - the quantitative accuracy stems from a 
fortuitous conspiracy of small numerical prefactors pq] . 



A. Berryon propagator 

Ultimately, we are interested in the properties of phys- 
ical electrons. To describe those we need to understand 
the properties of the electron-electron interaction medi- 
ated by the gauge field a^. To this end we proceed to 
calculate the Berry gauge field propagator by integrating 
out the fermion degrees of freedom from the Lagrangian 
(p8|). To one loop order this corresponds to evaluating 
the vacuum polarization bubble Fig. |^(a). Employing 
the standard rules for Feynman diagrams in the momen- 
tum space 1 39 the vacuum polarization reads: 



^Tr[Go(fc)7pGo(fc 



+ 9)7.]- (59) 
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Here G'o(fc) = ka'ja/k'^ is the free Dirac Green function, 
k = (fco,k) denotes the Euchdean three- momentum, and 
the trace is performed over the 7 matrices. 

The integral in Eq. ( ^9|) is a standard one (sec Ap- 
pendix for the details of computation) and the result is 



(60) 



where |g| = 1/(7^. The one loop effective action for the 
Berry gauge field therefore becomes (27r)~'^ / (fiqCs with 



N . 



2e2 



- ) af,{q)a^{-q) 



(61) 



At low energies and long wavelengths, \q\e'^ <C N/A, the 
fermion polarization completely overwhelms the original 
Maxwell bare action term and the Berryon properties be- 
come universal. In particular, the coupling constant 1/e^ 
drops out at low energies and only reappears as the ultra- 
violet cutoff. Physically, the medium of massless Dirac 
fermions screens the long range interactions mediated by 
a^. In QED3 this screening is incomplete: the gauge field 
becomes stiffer by one power of q but still remains mass- 
less, in accordance with our general expectations. This 
anomalous stiffness of justifies the quadratic level ex- 
pansion of Co ( ^sl) and renders higher order terms irrel- 
evant in the RG sense. The theory therefore clears an 
important self-consistency check. 

At low energies the fully dressed Berryon propagator 
is given as an inverse of the polarization. 



D^,iq) = U-^{q) 



(62) 



In order to perform this inversion we have to fix the 
gauge. To this end we implement the usual gauge fix- 
ing procedure by replacing q^q^/q'^ ^ (1 - C"^)?^?''/'?^ 
in Eq. (|60|). ^ > parametrizes the orbit of all covari- 
ant gauges. For example, ^ = corresponds to Lorentz 
gauge k^a^{k) = while ^ = 1 corresponds to Feynman 
gauge. Upon inversion we obtain the low energy Berryon 
propagator 



D^u[q) 



<5„.-^(l-0), (63) 



\q\N V q' 



in agreement with previous authors 



B. TF self energy and propagator 

TF propagator is a gauge dependent entity and one 
could therefore immediately object that as such it has 
no direct physical content and is of no interest. Such 
viewpoint, while expressed frequently, is actually quite 
naive. The reality is that in gauge theories various gauge- 
variant objects can often be connected to physical gauge- 
invariant quantities when computed within a particular 



choice of gauge. In practice, a rather typical occurrence 
is that a gauge-invariant physical propagator is given by 
a hugely non-local form which is basically impossible to 
compute except in a judiciously chosen gauge where it is 
related to a much simpler, and therefore far easier to com- 
pute, gauge-variant propagator. Consequently, a gauge- 
variant propagator computed along a particular gauge 
orbit often contains relevant information about the true 
dynamics of a gauge theory - the trick is to know how to 
extract this information. 

This general statement holds in the case of QED3 as 
well. TF propagator evaluated in an arbitrary covari- 
ant gauge parametrized by ^ contains useful information 
about the nature of the fermionic excitations of the sys- 
tem. We will show that its coupling to the massless 
gauge field destroys the usual Fermi liquid pole and re- 
sults in the propagator displaying a Luttinger-like be- 
havior, characterized by a small anomalous dimension. 
In the next subsection we argue that the physical elec- 
tron propagator of our theory, which can be related to a 
particular gauge invariant fermion propagator of QED3, 
exhibits the same Luttinger-like behavior. 

The lowest order self-energy diagram is depicted in Fig. 
|(b) and reads 



i 



E(fc) 



d^q 
(27r)3 



Df_,y{q)-ff_,Go{k + q)-fu 



(64) 



Again, the computation is rather straightforward (see 
Appendix for details) and the most divergent part is 



s(fc) 



4(2 ~ 30 
37r2iV 



/^In 



(65) 



where we have introduced the Feynman "slash" notation, 
= kf^jf^. 

To the leading 1/iV order, the inverse TF propagator 
is given by 



1 



with 



4(2 - 30 
37r2iV 



(66) 



(67) 



Higher order contributions in 1 /N will necessarily affect 
this result. Renormalization group arguments and 
non-perturbative approaches |43|| strongly suggest that 
Eq. (^6|) represents the start of a perturbative series that 
eventually resums into a power law: 



This implies real-space propagator of the form 



G(r) = A-"- 



•3+r; ■ 



(68) 



(69) 



16 



Thus, the TF propagator exhibits a Luttinger-like alge- 
braic singularity at small momenta, characterized by an 
anomalous exponent rj. In the Lorentz gauge = 0) 
we find ?7 = -8/3tt^N ~ -0.13, for N = 2. This 
rather small numerical value for the anomalous dimen- 
sion exponent (which is even considerably smaller for 
TV = 4 or = 6) indicates that the unraveling of the 
Fermi liquid pole in the original TF propagator brought 
about by its interaction with the massless Berry gauge 
field is in a certain sense "weak". Note also that rj is 
negative in the Lorentz gauge while it becomes positive, 
77 = 4/3tt'^N ~ 0.06 for N — 2, in the Feynman gauge 
= 1). The above results provide a strong indica- 
tion that the physical, gauge-invariant fermion propa- 
gator also has a Luttinger-like form, characterized by a 
small and positive anomalous dimension JTsf . We now 
show that this indeed is the case. 



C. Physical electron propagator 

Various spectroscopies on cuprates, such as ARPES 
and STM, as well as numerous optical and microwave 
techniques, all measure the spectral function of a real 
physical electron, not of TF. Therefore, we are ultimately 
interested in computing the propagator of the physical 
electron in our theory, 



(70) 



where ^(a;) is the original electron field operator appear- 
ing in (|l|) (the reader should recall that this operator 
already contains high-energy renormalizations built in at 
the very beginning). If we define a matrix 



n{x) 











(71) 



we can write G"^'°'^ in terms of TF fields as 

G°^''\x-x')^{h{x)^{x)^\x')h\x')) . (72) 

Ordinary spectroscopies reflect the diagonal part of the 
electron propagator 

[G"'"'=(a;-2;')]ii = {e-'^^^''''^-f^^'''^^{x)¥ {x')]i{) , 

(73) 

and similar expression for [G'^^'"^{x — a;')] 22- We may re- 
cast this in a more convenient form by writing the phase 
difference in the exponent as a line integral of a gradient 
along a straight line connecting x and x' , 



(74) 



and then expressing the phase gradient in terms of the 
two gauge fields and v^: 



[G<=^°^(x-x')]ii = (^e'^^'^^^+''^^'^'^'[^{x)¥{x')]ny 



(75) 



This expression only involves the coarse-grained Doppler 
and Berry gauge fields and TF, which are precisely the 
fields that enter our effective low-energy theory, and is 
therefore amenable to analysis. Note that (75) is only a 
long-distance approximation to the exact electron propa- 
gator ( [73I) which is defined through discrete vortex vari- 
ables entering via ipA (^b)- 

As discussed earlier, the Doppler gauge field is mas- 
sive in both normal and superconducting phases and 
therefore its fiuctuations will not affect the low energy, 
long wavelength properties of the electron propagator. 
We may thus remove it from the line integral Eq. ( |75| ) 
and focus on the quantity 



g{x ~ x') 



'*"T(a;)T(a;') 



(76) 



By considering the transformation properties of T(a;) un- 
der the gauge transformations with respect to it is 
easy to verify that Q(x — x') is gauge invariant. This 
quantity therefore represents a gauge invariant propa- 
gator of QED3 theory ( ^2[) and its knowledge allows us 
to reconstruct the diagonal components of the electron 
propagator by means of 



[G^'-(x-a;')]». = [i"foQ{x-x')]u 



(77) 



It seems natural to attempt to relate the components 
of the above gauge invariant propagator to the physical 
electron propagator, since the latter by definition must be 
gauge invariant under the transformations of the internal 
gauge field a^. 

The question arises how to evaluate the gauge invariant 
propagator Q(x — x'). This turns out to be a non-trivial 
issue since, despite its pleasing manifest invariance un- 
der gauge and spacetime symmetries, (^6|) also exhibits 
a severe (linear) ultraviolet divergence arising from the 
straight line integral of the gauge field. This renders it 
ill-defined in absence of some proper regularization (see 
below for more details). Here we adopt the approach dis- 
cussed by Brown |Q, ^ in the context of QED4. Brown 
proves that the following relation exists between gauge 
invariant propagator Q and gauge dependent propagator 
G H (see also Appendix C): 



^(r) =e^MG(r), 



(78) 



where r = x — x' and 



F{r) ^1 I d^z I <Pz'J^{z)D,,,{z - z')J,{z') , (79) 



with 



daS{z — x' — ar) 



(80) 



representing the source term for the line integral in Eq. 
([76|). In the above G(r) and D^^{r) refer to the real 
space gauge dependent fermion and Berryon propagators 
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respectively, obtained by Fourier transforming the ex- 
pressions (p^ ) and (|63|). By definition, both are to be 
computed in covariant gauge. Brown's result Eq. ( |7^ ) 
is an explicit statement of the fact that one can con- 
struct two gauge invariant propagators by using the line 
integral of the gauge field, Q and Q. This is a rather gen- 
eral feature of Abelian gauge theories [|5| and is easily 
generalized to QED3. Q and G can be formally related 
through a gauge transformation (see Appendix C) and 
one might think of Q as representing suitably regular- 
ized Q |47|. Alternatively, we can simply think of Q as 
being another QED3 fermion propagator invariant under 
gauge and spacetimc symmetries just like Q. We discuss 
the details pertaining to Eq. ( [zsl) in the Appendix C. 



To calculate g{r) from Eq. (78) we need to evaluate 
F{r). We proceed by first performing the z, z' integrals 
to obtain 



F{r) f da f dpr^D^, 
^ Jo Jo 



{r{a - I3))ry 



(81) 



By power counting the expression for F{r) suffers from 
linear UV divergence reflecting the singular behavior of 
the gauge field line integral at short distances. This sin- 
gularity is the main reason why direct computation of 
( [76| ) is such a frustrating task, despite its deceivingly 
compact and elegant form. A typical scheme to regularize 
this linear UV divergence interferes with gauge invariance 
and corrupts the effort to extract the true physical part 
of (|7^) which we expect to be scale-invariant. The advan- 
tage of Brown's approach ( |78| ) is twofold: it permits com- 
putation of a gauge invariant propagator Q in the covari- 
ant gauge where the UV divergence can be treated with 
the help of dimensional regularization, which respects the 
gauge invariance and preserves the long wavelength, low 
energy properties of the physical propagators, and the of- 
fending linear part of the UV divergence cancels out be- 
tween numerator and denominator in a gauge invariant 
manner (Appendix C). This allows us to extract a mean- 
ingful power law behavior for Q. To take advantage of di- 
mensional regularization we express V{r) = rf_iDf^u[r)r,j 
as a Fourier transform in d-dimensions. 



V{r) 



d'^k e 
N J (27r)° 



ik-r 



(1-0 



(fc-r)2 



fc2 



(82) 



and treat d as a continuous variable. This Fourier trans- 
form is evaluated in the Appendix, giving the result 



Vir) 



^^[i + id-2){i-oy-'. 



N 



(83) 



Substituting V{r{a — 0)) into Eq. (|8l|) and performing 
the remaining integrals by means of 



(84) 



we find, near d = 3 
F{r) = - 



4(2-^) hm 



Nn^ V 3 

4(2-0 



In(Ar) 



3-d 



(85) 



The UV divergence is now parametrized by the (r- 
independent) second term in the angular brackets. The 
leading long-distance behavior is contained in the log, 
implying a power law contribution e^''"^ cx r"**^^"^^/''^'^ 
to the gauge invariant propagator Q (r) . Combining ( |85| ) 
with (|9|) we obtain Q 



Sir) = A- 
or in the momentum space 



r3+V' 



(86) 



(87) 



where the anomalous dimension exponent rj' is given by 



4(2-0 
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37r27V 



(88) 



The last equation informs us that in the exponent rj' the 
gauge fixing parameter ^ has canceled out and Q{r) is 
indeed gauge invariant. We have thus verified, by explicit 
calculation to leading order in 1/7V, that Eq. (|7^ ) yields 
a gauge invariant TF propagator which we can connect 
to the physical electron propagator by means of Eq. (|7^ ) 

with g g. 

An interesting feature of the above result is that for 
^ = 2 we have F{r) = 0. In the QED literature this is 
known as "Yennie's gauge" and its significance is that in 
this particular gauge the diagonal components of the TF 
propagator are directly equal to g{r). In Yennie's gauge 
one can therefore evaluate various electron observables 
in terms of TF propagator without worrying about the 
exponential factors. This is just the situation we have 
anticipated in the previous subsection. One could also 
define an "anti-Yennie's gauge" (or a "non-local gauge" 
as it is known in the QED3 literature), ^ = |, in which 

7] = and the effect of the gauge field fluctuations on g{r) 
is contained entirely in F(r). To leading order in 1/A'^ 
this observation further solidifies the expectation that the 
leading log in the self energy indeed exponentiates and 
the low energy, long lengthscale propagator behaves as a 
power law. 

Another important feature to observe is that ry' > 
- the electron has acquired a positive anomalous dimen- 
sion. The positivity of rj' is mandatory from general con- 
siderations - once we perform the Euclidean rotation and 
obtain the real time electron propagator the conditions of 
unitarity and causality of our original problem demand 
77' > 0. This is also a physically sensible result implying 
that the interacting electron propagator (86) decays on 
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long lengthscales faster than the free BdG electron prop- 
agator. Interaction mediated by the massless gauge field 
destabilizes the Fermi liquid pole in the original propaga- 
tor and leads to Luttinger liquid-like power law fermionic 
correlator in the low energy, long wavelength limit - this 
is our algebraic Fermi liquid (AFL), a non- Fermi liquid 
symmetric phase which governs the physics of the pseu- 
dogap state. The positivity of 77' means that the interact- 
ing AFL propagator is less coherent than the free BdG 
electron propagator which is just what one expects on 
intuitive grounds. We therefore propose that G be iden- 
tified as the true electron propagator in our theory. 

The Luttinger-like electron propagator that follows 
from ds?! ) leads to a characteristic asymmetry between 
the energy and momentum distribution curves (EDC and 
MDC) observed in ARPES experiments ^ (note that, 
within our theory, this behavior is limited to the pseu- 
dogap state) - the MDC is a very sharp Lorentzian close 
to the Fermi surface while EDC is broad, reflecting de- 
coherence of physical electrons in the pseudogap state 
1^ ). This decoherence is relatively weak with the ex- 
plicit value of 77' ~ 0.27 for the case of an individual Cu02 
layer where N — 2. In bi-layer or multi-layer systems N 
could be four, six or even higher and 77' is even smaller. 
The reason for this increase in N is the anisotropy of the 
tunneling matrix element between the constituent Cu02 
planes within a multi-layer unit cell. This matrix ele- 
ment effectively vanishes near the nodes, along (±7r, ±7r) 
directions, but is appreciable elsewhere in the Brillouin 
zone . The result is that low energy BdG fermions on 
different constituent Cu02 layers remain decoupled while 
the vortex excitations on these same layers are strongly 
coupled within the unit cell of a multi-layer, since their 
coupling comes from an integral over the full Brillouin 
zone. This translates to a larger effective N in our QED3 
( ^ ) and to a corresponding reduction in the anomalous 
dimension: rj' ~ 0.13 for TV = 4 (bi-layer like YBCO) 
or r?' ~ 0.09 for iV = 6 (tri-layer like HgBazCaaCuaOs). 
The actual numerical value of the anomalous dimension 
exponent rj' is the "fingerprint" , a unique mathematical 
signature of the symmetric phase of QED3 and there- 
fore of the algebraic Fermi liquid (AFL) state within the 
pseudogap regime of underdoped cuprates. Determin- 
ing T]' directly from experiments, either through various 
spectroscopies or transport measurements, would be a 
major step toward testing theoretical ideas expressed in 
this paper. 

The exponent rj' of the physical electron propagator 
had come under much scrutiny as of late since several 
effective theories related to QED3 emerged recently in 
condensed matter physics, in problems like Heisenberg 
antiferromagnets or spin liquids. While in each case the 
physical content of these multiple reincarnations of QED3 
differs completely from the one discussed in the present 
paper and from each other, the issue of the gauge invari- 
ant QED3 fermion propagator and the value of rj' loom 
large in all these different contexts, for obvious reasons. 
In particular, 77' has been calculated recently by Rantner 



and Wen |5^ and also by Khveshchenko . The for- 
mer authors obtain 77' = -32/3Tr'^N || by performing a 
calculation of Q in the so-called axial gaug e [p5[ , in which 
the line integral of the gauge field in Eq. (|76D is taken to 
vanish for a particular direction in the real space. Neg- 
ative anomalous dimension 77' < would imply that the 
interacting electron propagator is more coherent at long 
distances than the propagator of a free electron and this 
is prohibited on general grounds, as discussed above. For 
example, negative 77' produces divergent electronic den- 
sity of states and leads to unphysical singular behavior in 
various thermodynamic and transport quantities. Thus, 
negative anomalous dimension for the physical electron 
should be, in our view, rejected out of hand. For reader's 
benefit, we should stress that we believe that the cal- 
culations carried out in Refs. ^ are perfectly 
correct, in the sense that rj' — -~32/3tt'^N indeed follows 
from the algebra once we adopt the axi al g auge regular- 
ization of ( [76| ) as implemented in ||5^, [52| and perform 
the calculation with logarithmic accuracy; we have done 
such a calculation ourselves and have obtained the same 
result. Operationally, the problem with computation of 
G (76) is not in the algebra but resides in the physical 
interpretation of the obtained results. In the axial gauge 
regularization followed by the momentum s pac e compu- 
tation employed by the authors of Refs. the 
negative value for 77' arises from the treatment of spuri- 
ous gauge singularities that are invariably introduced by 
writing the gauge boson propagator in the axial gauge 
p5| . These singularities then must be regularized in some 
way and this is done using an ad hoc principal value 
prescription for momentum space integrals. The prob- 
lem with this prescription is that the momentum space 
propagator in the axial gauge does not exist [ ^6| . Axial 
gauge is an example of a singularly non-covariant gauge 
(like Coulomb, temporal or similar gauges) and as such 
does not fully fix the gauge. The gauge transformations 
which are independent of the spacetime variable xi sin- 
gled out by the axial gauge but have arbitrary depen- 
dence on the remaining d — I variables (2:2, X3, . . .) are 
still allowed and cost no energy. Consequently, the gauge- 
variant two-point fermion propagator G(x, x') must van- 
ish whenever X2 X2, X3 x'^, . . . |Q. The same prob- 
lem resurfaces in a different form when one computes 
( [76| ) directly in the covariant gauge ||5^ : the expectation 
value of the transverse part (^ = 0) of the line integral, 

(exp(i Ofj^dSfj,)), decouples from the rest of the expres- 
sion. In fact, the result rj' = — 32/37r^A^ is easily under- 
stood as arising from the TF propagator in the Lorentz 
gauge G^=Q [rjL — — 8/37r^A^) being made more coherent 
through simply being multiplied by the said expectation 
value of the transverse part of line integral (77' = rjL + rjt 
where rjt = —8/tt'^N from the above computation of F, 

with ^ = 0). Since (exp(i a^(is^))^=o is an expecta- 
tion value of a phase factor, this is clearly a troubling re- 
sult - such a multiplicative factor can make the full prop- 
agator only less coherent than G^=o as |a; — a;'| — > 00. The 
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problem is that the transverse part of the hne integral is 
more divergent than just a simple logr appearing in the 
exponent. When the dominant (linear) divergence is in- 
cluded, the full propagator Q is exponentially suppressed 
(j] —^ +00) and cannot be computed without some phys- 
ically motivated UV regularization scheme ||4^ . 

We have also carried out our own calculation in the 
axial gauge using a different regularization scheme and 
obtained different (negative) exponent ?/ — — 16/37r^iV. 
We are thus forced to conclude that the axial gauge calcu- 
lations yield values of rj' that are regularization scheme 
dependent, and are therefore inherently unreliable. By 
contrast no such ambiguities arise when employing Eq. 
( [ZS] ) and the computation of Q yields physically reason- 
able positive anomalous dimension given by Eq. (|88|). 

We now discuss the off-diagonal (anomalous) compo- 
nents of the electron propagator. According to Eq. (|7l| ) 
we have 
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and similar expression for [G° (a; — a;')]2i. It is now less 
straightforward to interpret the phase factor in terms of 
our gauge fields and . One way to do this is to add 
and subtract ipAix') and write 

ipA{x) + ipb{x') = [ipAix) - ipA{x')] + [(Pa{x') + ipb{x')] 
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The first term on the right hand side is just like the one 
we encountered in our discussion of the diagonal prop- 
agator. The second term, however, involves a line in- 
tegral of Vf^ originating at an arbitrary fixed reference 
point in space-time. In the non-superconducting phase 
fluctuations in will clearly drive any such term to zero, 
making the off-diagonal terms of the electron propagator 
vanish. This is consistent with our general expectation 
that electron propagator does not exhibit anomalous off- 
diagonal correlations in the normal state. 

Finally, pulling different strands together, we can write 
down the full electron Green function in the AFL phase 
of the pseudogap state. To make connection with the 
notation prevalent in the condensed matter physics we 
perform a Euclidean rotation fcg — > ico in Eq. ( |87| ) and 
with help of Eq. (^^ we obtain 



G°'"'=(k,cj) = A- 



[e2 +A2 -c^2]l-„'/2' 



(91) 



where we have also restored the full electron dispersion 
Ek and the gap function Ak with the understanding that 
the above form for the propagator is strictly valid only in 
the vicinity of the nodal point and close to the isotropic 
limit. We note that the Eq. (91) implies an anomalous 
electron density of states 



(92) 



at low energies. It would be very interesting if such 
anomalous electron density of states could be measured 
in tunneling experiments. Simil arly , the Luttinger-like 
behavior of the propagator Eq. ( |9l| ) will be reflected in 
other physical observables. 



V. EFFECTS OF DIRAC ANISOTROPY IN 
SYMMETRIC QED3 

It is natural to examine to what extend is the theory 
modified by the inclusion of the Dirac anisotropy, i.e. 
the finite difference in the Fermi velocity vp and the gap 
velocity va ■ In the actual materials the Dirac anisotropy 
ao — ^ decreases with decreasing doping from ~ 15 in 
the optimally doped to '--^ 3 in the heavily underdoped 
samples. 

There are two key issues: first, for a large enough num- 
ber of Dirac fermion species N, how is the chirally sym- 
metric infrared (IR) fixed point modified by the fact that 
an 1, and second, as we decrease N, does the chiral 
symmetry breaking occur at the same value of N as in 
the isotropic theory. In this section we address in detail 
the first issue and defer the discussion of the second one 
to Part II. 

We determine the effect of the Dirac anisotropy, 
marginal by power counting, by the perturbative renor- 
malization group (RG) to first order in the large N ex- 
pansion. To the leading order S in the small anisotropy 
ao = 1 + S, we obtain the analytic value of the RG Pa^- 
function and find that it is proportional to S, i.e. in the 
infra red the a d decreases when 5 > and the anisotropic 
theory flows to the isotropic fixed point. On the other 
hand, when (5 < 0, ajj increases in the IR and again 
the theory flows into the isotropic fixed point. These re- 
sults hold even when anisotropy is not small as shown 
by numerical evaluation of the /3-function. Therefore, we 
conclude that the isotropic fixed point is stable against 
small anisotropy. 

Furthermore, we show that in any covariant gauge 
renormalization of S due to the unphysical longitudinal 
degrees of freedom is exactly the same along any space- 
time direction. Therefore the only contribution to the 
RG flow of anisotropy comes from the physical degrees 
of freedom and our results for (3ao stated above are in 
fact gauge invariant. 



A. Anisotropic QEDa 

In the realm of condensed matter physics there is no 
Lorentz symmetry to safeguard the space-time isotropy 
of the theory. Rather, the intrinsic Dirac anisotropy is 
always present since it ultimately arises from compli- 
cated microscopic interactions in the solid which even- 
tually renormalize to band and pairing amplitude dis- 
persion. Thus there is nothing to protect the difference 
in the Fermi velocity vp = and the gap velocity 
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I'A = from vanishing and, in fact, all HTS mate- 
rials are anisotropic. 

The value of a^i can be directly measured by the angle 
resolved photo-emission spectroscopy (ARPES), which is 
ultimately a "high" energy local probe of vp and va- 
Since QED3 is free on short distances, we can take the 
experimental values as the starting bare parameters of 
the field theory. 

The pairing amplitude of the HTS cuprates has dx^^y2 
symmetry, and consequently there are four nodal points 
on the Fermi surface with Dirac dispersion around which 
we can linearize the theory. Note that the roles of x and y 
directions are interchanged between adjacent nodes. As 
before, we combine the four two-component Dirac spinors 
for the opposite (time reversed) nodes into two four- 
components spinors and label them as (1,1) and (2,2) 
(see Fig. I). 

Thus, the two-point vertex function of the non- 
interacting theory for, say, 11 fermions is 

p(2)/ree ^ ^^^^ ^ Vp^lh + VAl2k2- (93) 

Therefore the corresponding non-interacting "nodal" 
Green functions are 

Here we introduced the (diagonal) "nodal" metric gj^l): 

(1) _ (2) _ , (1) _ (2) 2 (1) _ (2) _ 2 

900 — 900 — 9ll — 922 — ^F^ 922 — 9ii — 

well as the "nodal" 7 matrices 7". In what follows we 
assume that both vf and va are dimensionless and that 
eventually one of them can be chosen to be unity by an 
appropriate choice of the "speed of hght". 

Since 1 breaks Lorentz invariance of the theory 

and since it is the Lorentz invariance that protects the 
spacetime isotropy, we expect the /3 functions for ao to 
acquire finite values. However, the theory still respects 
time-reversal and parity and for N large enough the sys- 
tem is in the chirally symmetric phase. These symmetries 
force the fermion self-energy of the interacting theory to 
have the form 

^11 = A{hi, ^22) (7ofco + Wi?Ci7ifci + VAC2l2k2) ■ (95) 

where kii = k^g^^Jk^ and k22 = k^^g'l^}ki,. The coeffi- 
cients Ci are in general different from unity. Furthermore, 
there is a discrete symmetry which relates flavors 1, 1 and 
2, 2 and the x and y directions in such a way that 

S22 = ^(^22, fcll) (7ofco + (271^1 + VFC,ll2k2) ■ (96) 

In the computation of the fermion self-energy, this dis- 
crete symmetry allows us to concentrate on a particular 
pair of nodes without any loss of generality. 

B. Gauge field propagator 

As discussed above in the isotropic case, the effect of 
vortex-antivortex fluctuations at T=0 on the fermions 



can, at large distances, be included by coupling the nodal 
fermions minimally to a fluctuating U{1) gauge fleld with 
a standard Maxwell action. Upon integrating out the 
fermions, the gauge field acquires a stiffness proportional 
to fc, which is another way of saying that at the charged, 
chirally symmetric fixed point the gauge field has an 
anomalous dimension tja — ^ (for discussion of this point 
in the bosonic QED see p8| ). 

We first proceed in the transverse gauge (A;^a^ = 0) 
which is in some sense the most physical one considering 
that the V x a is physically related to the vorticity, i.e. 
an intrinsically transverse quantity. We later extend our 
results to a general covariant gauge. To one-loop order 
the screening effects of the fermions on the gauge field 
are given by the polarization function 

nM.(^) = y E J J^Tr[GUqKGlS{q + k)^:] 

(97) 

where the index n denotes the fermion "nodal" flavor. 
The above expression can be evaluated straightforwardly 
by noting that it reduces to the isotropic H^i,(fc) once 
the integrals are properly rescaled ||l^. The result can 
be conveniently presented by taking advantage of the 
"nodal" metric as 

(98) 

Note that this expression is explicitly transverse, 
fc^H^y(fc) — H^^(fc)fcj^ — 0, and symmetric in its space- 
time indices. It also properly reduces to the isotropic 
expression when vp — va = 1- 

However, as opposed to the isotropic case, it is not 
quite as straightforward to determine the gauge field 
propagator For example, as it stands the polar- 

ization matrix ( p8| ) is not invertible, which makes it nec- 
essary to introduce some gauge-fixing conditions. In our 
case the direct inversion of the 3x3 matrix would ob- 
scure the analysis and therefore, we choose to follow a 
more physical and notationally transparent line of rea- 
soning which eventually leads to the correct expression 
for the gauge field propagator. Upon integrating out the 
fermions and expanding the effective action to the one- 
loop order, we find that 

Ceff[a^] = (H^S + H^.)a^a, (99) 

where the bare gauge field stiffness is 

- (^.^ - ^) • (100) 

At this point we introduce the dual field 5^ which is re- 
lated to Ufj, as 

bf_t = efi^xqi^ax- (101) 
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Physically, the 6^ field represents vorticity and we inte- 
grate over all possible vorticity configurations with the 



restriction that 6^ is transverse. We note that 



(102) 



where X/^'s are functions of and upon a straightforward 
calculation they can be found to read 
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N 



2e^ WvpVA 
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n=l,2 



(103) 



where ^ ^ u X E {0,1,2}. At low energies we 
can neglect the non-divergent bare stiffness and thus set 
1/e^ = in the above expression. 

The expression (102) is manifestly gauge invariant and 
has the merit of not only being quadratic but also diago- 
nal in the individual components of 5^. Thus, integration 
over the vorticity (even with the restriction of transverse 
b^) is simple and we can easily determine the b^ field 
correlation function 



(bA) = — 




(104) 



The repeated indices are not summed in the above ex- 
pression. Note that, in addition to being transverse, 
{bfjbu) is also symmetric in its space time indices. 

It is now quite simple to determine the correlation 
function for the field and in the transverse gauge we 
obtain 



(bjbi). 



(105) 



Using the transverse character of {bfj^y) (which is inde- 
pendent of the gauge) the above expression can be further 
reduced to 

D,M = ^ - ^) (6^) - {b,b.)^ . (106) 

It can be easily checked that in the isotropic limit the 
expression (106) properly reduces to the results obtained 
in a different way. 

We can further extend this result to include a general 
gauge by writing 



1 



{b') - {bM 



(107) 



where ^ is our continuous parameterization of the gauge 
fixing. This expression can be justified by the Fadeev- 
Popov type of procedure starting from the Lagrangian 



/• — ( n 1 2fc^ fc^fc^ 



(108) 



where the stiffness for the unphysical modes was judi- 
ciously chosen to scale as A: in a particular combination 



of the physical scalars of the theory. Note that (6^) can be 
determined without ever considering gauge fixing terms. 
In this way, the extension of a transverse gauge ^ = 
to a general covariant gauge is accomplished by a sim- 
ple s ubstitution k^ky — > (1 — |)fc^fcy. The expression 
(107) is our final result for the gauge field propagator in 
a covariant gauge. 



C. TF self energy 

As discussed in Sec. IV, S is not gauge invariant in 
that it has an explicit dependence on the gauge fixing 
parameter. As we will show in this Section (and more 
generally in the Appendix D), the renormalization of S 
by the unphysical longitudinal degrees of freedom does 
not depend on the space-time direction: the term in S 
which is proportional to 70 is renormalizcd the same way 
by the gauge dependent part of the action as the terms 
proportional to 71 and 72. Therefore, the only contri- 
bution to the RG flow of ajj comes from the physical 
degrees of freedom. 

We denote the topological fermion self-energy at the 
node n by Tjn{q)- Hence, to the leading order in large N 
expansion we have 



S„(q) : 
or explicitly 
S„(q) = 



(Pk 
(2^ 



l';G^i{q-k)^:D^,(k). (109) 



(fk 



(9-fc)A7r 



(2^)3 "'(<z~%5^„(9-fc), 



-7"^M^(fc), 



(110) 



where the gauge field p ropa gator D^^, is already screened 
by the nodal fermions (107). Using the fact that 



It^lxlv = i^tiXulbli + 



where /x, i', A € {0,1,2} and 75 
easily see that 



'iixl,^ - 5^,^"f\ + 5xul^l, (111) 
-170717273, we can 



(112) 



where we used the symmetry of the gauge field propaga- 
tor tensor D^^. Thus, 



^n{q) 



d'k {q-k)x{2glp:-^^g-^,)D^,{k) 



iq-k)f,gjl^{q-k)^ 



(113) 



and as shown in the Appendix D at low energies this can 
be written as 



S„(g)=-Er;;(7;9M)ln 



A 



y^<lc.9apQP 



(114) 



Here A is an upper cutoff and the coefficients ij are func- 
tions of the bare anisotropy, which have been reduced 
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to a quadrature (see Appendix D). It is straightforward, 
even if somewhat tedious, to show that in case of weak 
anisotropy (uf = 1 + ua = 1), to order , 

al 8 3 1 2 



1- 2^-^(40 -70 (115) 
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In the isotropic hmit [vp = wa = 1) we regain rf^ 
—8(1 — I^/Stt^A^ as previously found by others. 



D. Dirac anisotropy and its /3 function 

Before plunging into any formal analysis, we wish to 
discuss some immediate observations regardin g th e RG 
flow of the anisotropy. Examining the Eq. (114) it is 
clear that if 77" = r]^ then the anisotropy does not flow 
and remains equal to its bare value. That would mean 
that anisotropy is marginal and the theory flows into the 
anisotropic fixed point. In fact, such a theory would have 
a critical line of ao- For this to happen, however, there 
would have to be a symmetry which would protect the 
equality 77" = ■ For example, in the isotropic QED3 
the symmetry which protects the equality of ry's is the 
Lorentz invariance. In the case at hand, this symmetry 
is broken and therefore we expect that 77" will be differ- 
ent from 772 , suggesting that the anisotropy flows away 
from its bare value. If we start with a/j > 1 and find that 
?72^ > at some scale p < A, we would conclude that 
the anisotropy is marginally irrelevant and decreases to- 
wards 1. On the other hand if 772^ < 'q\^ , then anisotropy 
continues increasing beyond its bare value and the theory 
flows into a critical point with (in)finitc anisotropy. 

The issue is further complicated by the fact that 77^ 
is not a gauge invariant quantity, i.e. it depends on 
the gauge fixing parameter ^. The statement that, say 
> makes sense only if the ^ dependence of 77" 
and 772 is exactly the same, otherwise we could choose 
a gauge in which the difference 772 — tfi can have either 



sign. However, we see from the equations ( 115 -117) that 
in fact the ^ dependence of all 77's is indeed the same. Al- 
though it was explicitly demonstrated only to the 0{S'^), 
in the Appendix D we show that it is in fact true to 
all orders of anisotropy for any choice of covariant gauge 
fixing. This fact provides the justification for our proce- 
dure. Now we supply the formal analysis reflecting the 
above discussion. 

The renormalized 2-point vertex function is related to 
the "bare" 2-point vertex function via a fermion field 
rescaling as 



(118) 



It is natural to demand that for example at nodes 1 and 
1 at some renormalization scale p, Tj^ (p) will have the 
form 



^R^P) = 70P0 + Vp^ipi + Ua72P2- 



(119) 



Thus, the equation (119) corresponds to our renormaliza- 



tion condition through which we can eliminate the cutoff 
dependence and calculate the RG flows. 
To the order of 1 /N we can write 



(120) 



where we used the fermionic self-energy ( |114| ) . Multiply- 
ing both sides by 70 and tracing the resulting expression 
determines the field strength renormalization 
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(121) 



We can now determine the renormalized Fermi and gap 
velocities 



VF 



and 



VA 



(1 - 77II In ^)(1 + 771" In ^) « 1 - (77^^ - 771") In ^ 

(122) 



(l-77iMn-)(l-K77iMn-) 
P P 



l-(^o -%'')ln 



A 



P 

(123) 



The corresponding renormalized Dirac anisotropy is 
therefore 



^«a.,(l-(772"-7;,")ln^) 



The RG beta function can now be determined 



" d\np 



(124) 



(125) 



In the case of weak anisotropy {vp = 1 + 5, wa = 1) 
the above expression can be determined analytically as 
an expansion m S. Using Eqs.( pT6HTT^ ) we obtain 



3tt^N V 5 



(126) 



Note that this expression is independent of the gauge fix- 
ing parameter ^. For < 5 ^ 1 the /3 function is positive 
which means that anisotropy decreases in the IR and thus 
the anisotropic QED3 scales to an isotropic QED3. For 
— 1 <C (5 < the P function is negative and in this case 
the anisotropy increases towards the fixed point ao — 1, 
i.e. again towards the isotropic QED3. Note that for 
(5 > 2, /3 < which may naively indicate that there is a 
fixed point at 5 = 2; this however cannot be trusted as it 
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FIG. 4: The RG /3-function for the Dirac anisotropy in units 
of 8/3tv'^N. The solid line is t he numerical integration of 
the quadrature in the Eq. (D8) while the dash-dotted line 
is the analytic al expansion around the small anisotropy (see 
Eq. ( 115| - 117|) ). At ao ~ 1, Pao crosses zero with posi- 
tive slope, and therefore at large length-scales the anisotropic 
QED3 scales to an isotropic theory. 



is outside of the range of validity of the power expansion 
of 77^ . The numerical evaluation of the quadrature in Eq. 
( |D§| ) shows that, apart from the isotropic fixed point and 
the unstable fixed point at au = 0, /?q,q does not vanish 
(see Fig. |4|). This indicates that to the leading order in 
the 1/iV expansion, the theory flows into the isotropic 
fixed point. 



VI. SUMMARY AND CONCLUSIONS 

By appealing to the unique features of high- Tc cuprates 
- strong electron correlations, unconventional order pa- 
rameter symmetry and pronounced fluctuation effects - 
we argued in favor of the inverted approach to the prob- 
lem of describing various thermodynamic phases appear- 
ing in the underdoped region. This inverted approach 
can be thought of as a "Fermi liquid" theory of the phase 
fluctuating d-wave superconductor where the role of the 
Fermi energy as the large energy scale of the problem 
is played by A, the amplitude of the pseudogap, which 
we assume to be predominantly of the pairing nature. 
Under the umbrella of this pseudogap - inside the pair- 
ing protectorate - we identify the BdG quasiparticles as 
the relevant low lying fermionic excitations of the theory 
and study their evolution under the effects of interactions 
mediated by vortex-antivortex fluctuations. By carefully 
treating these interactions we find that the low energy 
effective theory for the quasiparticles inside the pairing 
protectorate is the (2+1) dimensional quantum electro- 
dynamics (QED3) with inherent spatial anisotropy, de- 
scribed by Lagrangian Co specified by Eqs. ( ^TUsT] ). 

Within the superconducting state the gauge fields of 
the theory are massive by virtue of vortex defects being 
bound into finite loops or vortex-antivortex pairs. Such 
massive gauge fields produce only short ranged interac- 
tions between our BdG quasiparticles and are therefore 
irrelevant: in the superconductor quasiparticles remain 
sharp in agreement with prevailing experimental data 



p9[ . Loss of the long range superconducting order is 
brought about by unbinding the topological defects - 
vortex loops or vortex-antivortex pairs - via Kosterlitz- 
Thouless type transition and its quantum cousin. Re- 
markably, this is accompanied by the Berry gauge field 
becoming massless. Such massless gauge field mediates 
long range interactions between the fermions and be- 
comes a relevant perturbation. Exactly what is the conse- 
quence of this relevant perturbation depends on the num- 
ber of fermion species N in the problem. For cuprates 
we argued that N — 2ncuO where ncuO is the number 
of CuOa layers per unit ceU. If < A^^ =i 3 the 
interactions cause spontaneous opening of a gap for the 
fermionic excitations at T = 0, via the mechanism of 
chiral symmetry breaking in QED3 |l9| . Formation of 
the gap corresponds to the onset of AF SDW instability 
1^ , p2| which must be considered as a progenitor of the 
Mott-Hubbard-Neel antiferromagnet at half-filling. If, on 
the other hand, N > Nc as will be the case in bilayer 
or trilayer materials, the theory remains in its chirally 
symmetric nonsuperconducting phase even as T — > 
and AF order arises from within such state only upon 
further underdoping (Fig. 1). We call this symmetric 
state of QED3 an algebraic Fermi liquid (AFL). In both 
cases AFL controls the low temperature, low energy be- 
havior of the pseudogap state and in this sense assumes 
the role played by Fermi liquid theory in conventional 
metals and superconductors. In AFL the quasiparticle 
pole is replaced by a branch cut - the quasiparticle is no 
longer sharp - and the gauge invariant electron propaga- 
tor acquires a Luttinger-like form Eq. (^7|) with positive 
anomalous dimension rj' = IB/Btt^TV. To our knowledge 
this is one of the very few cases where a non-Fermi liq- 
uid nature of the excitations has been demonstrated in 
dimension greater than one in the absence of disorder or 
magnetic field. This Luttinger-like behavior of AFL will 
manifest itself in anomalous power law functional form 
of many physical properties of the system. 

Dirac anisotropy, i.e. the fact that vp 7^ t^A, plays im- 
portant role in the cuprates where the ratio ao = vf/va 
in most materials ranges between 3 and 15—20. Such 
anisotropy is nontrivial as it cannot be rescaled and it 
significantly complicates any calculation within the the- 
ory. Using the perturbative renormalization group theory 
we have shown that anisotropic QED3 flows back into 
isotropic stable fixed point. This means that for weak 
anisotropy at long lengthscales the universal properties 
of the theory are identical to those of the simple isotropic 
case. It remains to be seen what are the properties of 
the theory at intermediate lengthscales when anisotropy 
is strong. 

Our theory of the pseudogap state gets its inspiration 
and builds on the ideas originally articulated by Emery 
and Kivelson in Ref. |^ and by Randeria and collabora- 
tors . These ideas were later explored and extended in 
various directions by others ||l|, |l^, |l|, |62[ g . These 
approaches share a common philosophical platform of as- 
suming that the pseudogap is primarily due to pairing in 
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the p-p channel and in the underdoped regime supercon- 
ducting long range coherence is destroyed by the phase 
fluctuations. The essential differences lie in the imple- 
mentation of these ideas. Balents et al. [Q, for in- 
stance, argue for true separation of spin and charge | |3^ 
in their " nodal liquid" phase which then gives way to an 
unconventional antiferromagnet "AF*" with deconfined 
spin 1/2 excitations. Our theory also shares consider- 
able formal similarities with the SU(2) gauge theories of 
Wen, Lee and co-workers |^ and works by Aitchinson 
and Mavromatos in that their low energy effective 
theory is also related to QED3. In those cases, however, 
the physical content of the theory - the identity of its 
excitations - are totally different. 

Finally there is a class of theories where the pseudogap 
is assumed to be due to some competing order, usually in 
the p-h channel. This includes SO (5) Q, d-density wave 
jSSt and various other competing orders |6^ . At present 
time we believe that experimental evidence favors the 
pairing origin of the pseudogap, however the evidence 
is far from conclusive and experiments can be found to 
support virtually any aspect of the above mentioned the- 
ories. We believe, therefore, that the road ahead necessi- 
tates making specific predictions based on controlled and 
well defined approaches. The QED3 theory of the pair- 
ing pseudogap, as presented in this paper, starts from 
a remarkably simple set of assumptions, and via manip- 
ulations that are controlled in the sense of expan- 
sion, arrives at nontrivial consequences, including the al- 
gebraic Fermi liquid and the antiferromagnet. Another 
notable feature of this theory is that it is fully calculable: 
having the explicit form of the low energy effective action, 
free of constraints and uncertainties, physical observables 
can be computed through a systematic procedure |67|. 



excitations relevant for low temperatures (T ^ T* 
in the underdoped regime. 



deep 



2D thermal vortex-antivortex fluctuations 



In order to perform specific computations we have to 
adopt a model for vortex-antivortex excitations. We 
will use a 2D Coulomb gas picture of vortex-antivortex 
plasma. In this model (anti)vortices are either point-like 
objects or are assumed to have a small hard-disk radius 

of size of the coherence length which emulates the core 

_i 

region. As long as ^0 ^ 2 , where n = riy + Ua is the 
average density of vortex defects, the two models lead to 
very similar results and both undergo a vortex-antivortex 
pair unbinding transition of the Kosterlitz-Thouless va- 
riety. 

Above the transition we have 



exp [-/3 / d\Co] = 2-^' I ^^(^) (Al) 

•' A,B '' 

x,5[V X V - iV X (VipA + Vv^s)] 
X(5[V X a- iy X (V(pA - V(^i3)] . 



The phase tp{r) is due solely to vortices and we can 
rewrite (Al) as: 
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APPENDIX A: JACOBIAN £0 

Here we derive the explicit form of the "Jacobian" Cq 
for two cases of interest: i) the thermal vortex-antivortex 
fluctuations in 2D layers and ii) the spacetime vortex loop 



E 



NJNJ 

N^,Na A,B i 

N„ 



Na 



c4p,(r) - E<5(r - r,)]<5[pa(r) -Y.^^^- r.)] 

i 3 

x5[h{r) - Trg 5(r - rf ) + tt^ <5(r - rf ) 

i i 



Here Ny{Na) is the number of free vortices (antivortices), 
Ni — Ny + Na, Ti (rj) are vortex (antivortex) coordi- 
nates and Pi,(r) (pa(r)) are the corresponding densities. 
6(r) = (V X a(r)), = nip^ - pB _ pA ^ pB) ^nd 
is the core energy which we have absorbed into Cq for 
convenience. We now express the above (5-functions as 
functional integrals over three new flelds: d„(r), c?a(r) 



25 



and K(r): 



where 



E 



2-N, 



NJNJ 

N^,Na A,B i 



J VdyVdaVnexpli J (frdy{py{r) - ^(5(r - r^)) 

i 

Na 

d\dMv)-Y,5{v~v,)) 



+i 







N„ 



N„ 



The integration over iJ-functions in the exponential is eas- 
ily performed and the summation over A(B) labels can 
be carried out explicitly to obtain: 



E 



VdyVda'DK exp 



i / d^r(dypy + daPa + Kb){AA) 



Y]^ / d^nexp{-idy{ri))cos{nK{rij) 



Y[ / rf^rj cxp(-irfa(i"j))cos(7rK(rj)). 



In the thermodynamic limit the sum (A4) is dominated 
by Nyi^a) = {Nv{a)), where (-/V„(q)) is the average num- 
ber of free (anti)vortices determined by solving the full 
problem. Furthermore, as (-/V„((j)) ^ oo in the thermody- 
namic limit the integration over dy{r), da(j) and K(r) can 
be performed in the saddle-point approximation leading 
to the following saddle-point equations: 

- py{r) + {Ny)ny{r) = (A5) 

~Pa{r) + {Na)na{r) = (A6) 

-bir) + [{Ny)ny{r) + (Na)nair)] (A7) 
XTT tanh(7rK(r)) = 



with 



Qd.„,{r) cosh(7rK(r)) 
/ d2^/e'im(r')cosh(7rK(r'))' 



m = a.v. 



Eqs. ( A5|- at]) follow from functional derivatives of (A4) 
with respect to d„(r), da{v) and ^(r), respectively. We 
have also built in the fact that the saddle-point solutions 
occur at dy idy, da ida, k in. 



The saddle-point equations ( A5 ■ A7 ) can be solved ex- 
actly leading to: 

dy{r) = lnpy{r) — lncosh(7rK(r)) (A8) 
da{r) = Inpa(r) — lncosh(7rK(r)) (A9) 



K{r) — — tanh 



6(r) 



n{py{r) + Pair)) 



(AlO) 



Inserting (|A8yAlO|) back into (0) finally gives the en- 
tropic part of Co/T: 



Py In py + Pain Pa (Vxa)^tanh 



(V X a) 



+ (Pv + Pa) In cosh tanh 



(V X a), 

T^iPv + Pa) 



T^iPv + Pa) 

(All) 



where p„(a)(r) are densities of free (anti)vortices. We 
display Cq in this form to make contact with familiar 
physics: the first two terms in (All) are the entropic 
contribution of free (anti)vortices and the Doppler gauge 
field Vxv ^ 7r(/9„— Pa) ((Vxv) =0). The last two terms 
encode the "Berry phase" physics of topological frustra- 
tion. Note that the "Berry" magnetic field 6 = (V x a)2 
couples directly only to the total density of vortex defects 
Pv + Pa and is insensitive to the vortex charge. This is 
a reflection of the Z2 symmetry of the original problem 
defined on discrete (i.e., not coarse-grained) vort ices. We 
write Pt,(a)(r) = {Pv{a)) + Spv{a){'r) and expand ( |AllD to 
leading order in 5pyi^a) and V x a: 

Co/T ^ (V X v)V(2^2^0 + (V X a)V(2^'n,), (A12) 

where n/ = {py) + {pa) is the average density of free vortex 
defects. Both v and a have a Maxwellian bare stiffness 
and are massless in the normal state. As one approaches 
7c, TLi ~ ^g^2 ~f 0, where i,sc{x, T) is the superconducting 
correlation length, and v and a become massive (see the 
main text). 



2. Quantum fluctuations of (2+l)D vortex loops 



The expression for /lobM] given by Eq. (|39|) follows 
directly once the system contains unbound vortex loops 
in its ground state and thus can respond to the exter- 
nal perturbation A^^^'^ over arbitrary large distances in 
(2+l)-dimensional spacetime. This is already clear at 
intuitive level if we just think of the geometry of infinite 
versus finite loops and the fact that only unbound loops 
allow vorticity fluctuations, described by {ju,{q)jv(~q)), 
to proceed unhindered. Still, it is useful to derive (|39|) 
and its consequences belabored in Section II within an 
explicit model for vortex loop fluctuations. 

Here we consider fluctuating vortex loops in continu- 
ous (2-1- 1)D spacetime and compute £o[juJ_using duality 
map to the relativistic Bose superfluid Q. We again 
start by using our model of a large gap s-wave supercon- 
ductor which enables us to neglect in the fermion ac- 
tion and integrate over it in the expression for C{)[v^, a^] 
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ll). This gives: 



oc ^ N 



m 

N=0 1=1 



where 



5 = ^ S'o / ds; + i ^ i dsi (j) dsi,g{\ 
1=1 Li'=i-' 



(A13) 

xi[si] - xi>[si-]\) 
(A14) 



and 



AT 



(A15) 

In the above equations N is the number of loops, si is 
the Schwinger proper time (or "proper length" ) of loop 

S'o is the action per unit length associated with motion 
of vortex cores in (2+l)-dimensional spacetime (in an 
analogous 3D model this would be Ec/T, where £c is the 
core line energy), g(|a;;[s;] — a;;/[s;']|) is the short range 
penalty for core overlap and L denotes a line integral. We 
kept our practice of including core terms independent of 
vorticity into Cq. Note that vortex loops must be periodic 
along T reflecting the original periodicity of (p(y,t). 

We can think of vortex loops as worldline trajectories 
of some relativistic charged (complex) bosons (charged 
since the loops have two orientations) in two spatial di- 
mensions. £o without the (5-function describes the vac- 
uum Lagrangian of such a theory, with vortex loops rep- 
resenting particle-antiparticle virtual creation and anni- 
hilation process. The duality map is based on the follow- 
ing relation between the Green function of free charged 
bosons and a gas of free oriented loops: 



G{x) = mwi^)) 



(27r)3 A;2 + 



dse 



(27r)3 



dse I I e 

V4^s; 



(A16) 



where nid is the mass of the complex dual field (l){x). 
Within the Feynman path integral representation we can 
write: 



(47rs) 



3/2 px(s)=x 
12/ I ^ ' 

— -X s 



x(0)=0 



'Dx{s') exp 



1 



ds'x'is') 



(A17) 

where x = dx/ds. Furthermore, by simple integration 



( A16| ) can be manipulated into 

^e-^™^^I?a:(s')exp ds'x^is') 



(A18) 

where the path integral now runs over closed loops 
(x(s) — x{Q)). In the dual theory this can be reexpressed 

as 

^ln{e+ml). (A19) 



Combining ( A18 ) and ( A19 ) and using 

Tr [ln(-a2 + mj)] = lnDet(-52 + ml) = W (A20) 
finally leads to: 



N r 



N=0 ' 1 = 1 



exp 



°° ds_ 

Si 
N 



1=1 •'^ 



(A21) 



This is nothing else but the partition function of the free 
loop gas. The size of loops is regulated by vrid- As rud — > 
the averag e loo p size diverges. On the other hand, 
through W ( |A2C| ), we can also think of ( |A2l[ ) as the 
partition function of the free bosonic theory. 
To exploit this equivalence further we write 



where 



Cd=m^ + m'M^ + ^-\cj,\' 



(A22) 



(A23) 



and argue that Zd describes vortex loops with short range 
interactions in ( [AT^ ). This can be easily demonstrated 
by decoupling through Hubbard-Stratonovich trans- 
formation and retracing the above steps. The reader is 
referred to the book by Kleinert for further details of the 

above duality mapping 

We can now rewrite the (5-function in ( A13| ) as 



.fc(.)-,v(.))-./l>»,exp(,/<i».„,0,-„J) 

(A24) 

and observe that in the above language of Feynman path 



integrals in proper time i J d'^XK^n^ (A15) assumes the 



meaning of a particle current three-vector coupled to 
a three- vector pot entia l k^. Employing the same argu- 
ments that led to ( A23| ) we now have: 



Cd[Ktj\ = 1(9 - in)(t)Y 



(A25) 
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which leads to 



(A26) 

In the pseudogap state vortex loop unbinding causes loss 
of superconducting phase coherence. In the dual lan- 
guage, the appearance of such infinite loops as to^ 
implies superflu idity in the system of charged bosons de- 
scribed by £d ( |A23|) . The dual off-diagonal long range 
order (ODLRO) in {(j){x)(j)* {x')) means that there are 
worldline trajectories that connect points x and x' over 
infinite spacetime distances - these infinite worldlines are 
nothing else but unbound vortex paths in this "vortex 
loop condensate" . Thus, the dual and the true super- 
conducting ODLRO are two opposite sides of the same 
coin. 

Cd[Kfj] is the Lagrangian of this dual supcrfluid in pres- 
ence of the external vector potential k^. In the ordered 
phase, the response of the system is just the dual version 
of the Meissner effect. Consequently, upon functional in- 
tegration over we are allowed to write: 



= / Vk^s' 



(A27) 

where = Kf/))]^, with (0) being the dual order pa- 
rameter. The remaining functional integration over 
finally results in: 



(A28) 



We have tacitly assumed that the system of loops is 
isotropic. The intrinsic anis otrop y of the (2+l)D the- 
ory is easily reinstated and (A28) becomes Eq. ( p9| ) of 
the main text. 

It is now time to recall that we are interested in a d- 
wave superconductor. This means we must restore 
to the problem. To accomplish this we engage in a bit 
of thievery: imagine now that it was whose coupling 
to fermions was negligible and we could integrate over it 
in (pT|). We would then be left with only the ^-function 
containing a^;. Actually, we can compute such £o[&/i/7r], 
where d x da = b, without any additional work. Note 
that Co contains only vorticity independent terms. We 
can equally well proclaim that it is the A{B) labels that 
determine the true orientation of our loops while the ac- 
tual vorticity is simply a gauge label - in essence, and 



no right to coarse-grain to begin with and would have 
to face up to its purely Z2 character (see Section II). Still, 
the above reasoning does illustrate that the Maxwellian 
stiffness of follows the same pattern as that of v^: 
both are determined by the order parameter {(f)) of con- 
densed dual superfluid. Thus, we can write the correct 
form of Cq, with both and fully included into our 
accounting, as 



(d X v)^{d X v)f, {d X a)f,{d X a)^ 



47r2|(0)| 



(A30) 



where our ignorance is now stored in computing the ac- 
tual value of (0) from the original parameters of the d- 
wave superconductor problem. With anisotropy restored 
this is precisely Eq. (|45|) of Section II. 



APPENDIX B: FEYNMAN INTEGRALS IN QED3 

Many of the integrals encountered here are consid- 
ered standard in particle physics. Since the techniques 
involved are not as common in the condensed matter 
physics we provide some of the technical details in this 
Appendix. A more in-depth discussion can be found in 
many field theory textbooks p9|. 



1. Vacuum polarization bubble 

The vacuum polarization Eq. ( ^ ) can be written more 
explicitly as 



U^,{q) = 27VTr[7 

a ^fi '^(3 



with 



cPk ka{kp + qp) 
(27r)3 k^{k + q)^ ' 



(Bl) 



(B2) 



The integrals of this type are most easily evaluated by 
employing the Feynman parametrization [ ^9[ . This con- 
sists in combining the denominators using the formula 



A-Bb r(a)r(6) io """^ 



b-1 



x^-^jl-x) 
[xA + {1- x)B]''+'> 



(B3) 



valid for any positive real numbers a, 5, A, B. Setting 
A — k'^ and B = {k + qY allows us to rewrite ( p^ ) as 



obtain: 



1 vorticity IS simply a gauge laoei - m essence, ana „i „ 

trade places. After the same algebra as before we lapiq) = dx - — 



kaikfj + qp) 



(A29) 



which is just the Maxwell action for a^. 

Of course, this simple argument that led to ( A29) is il- 
legal. We cannot just forget v^. If we did we would have 



)3 [{k + {l-x)q)^ +x{l~x)q^]^' 

(B4) 

We now shift the integration variable k ^ k — {1 — x)q 
to obtain 

, f d^k kakp - x{l - x)qaqp . . 
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where we have dropped terms odd in k which vanish by 
symmetry upon integration. We now notice that k^kp 
term will only contribute if a = /? and we can therefore 
replace it in Eq. (B5) by ^Sapk'^. With this replacement 
the angular integrals are trivial and we have 



Ia0{q) = [ dx [ dkk 
Jo Jo 



2 l'>apk^ - x{l - X)qaq(3 



[fc2 + x{l - x)q'^] 



,212 



(B6) 



The only remaining difficulty stems from the fact that the 
integral formally diverges at the upper bound. This di- 
vergence is an artifact of our linearization of the fcrmionic 
spectrum which breaks down for energies approaching the 
superconducting gap value. It is therefore completely 
legitimate to introduce an ultraviolet cutoff. Such UV 
cutoffs however tend to interfere with gauge invariance 
preservation of which is crucial in this computation. A 
more physical way of regularizing the integral Eq. (B6) 
is to recall that the gauge field must remain massless, i.e. 
ILn,{q ^ 0) = 0. To enforce this property we write Eq. 
as 

n,,,(<z) = 27VTr[7„7,,7/37.][/a/3(g) - 1^0)1 (B7) 



and we see that proper regularization of Eq. (B6) involves 
subtracting the value of the integral a,t q — 0. The re- 
maining integral is convergent and elementary; explicit 
evaluation gives 



Iap{q) 



(B8) 



Inserting this in (B7) and working out the trace using 
Eqs. (^) and (|5|) we find the result (|60|). Identical result 
can be obtained using dimensional regularization. 



Combining the denominators using Eq. (B3) and follow- 
ing the same steps as in the computation of polarization 
bubble above we obtain 

C{k) = \j dx^ J d^q (B12) 

{2x - l)(fc ■ q)^ - (1 - x)k'^q^ + xjl - xfk^ 
^ [g2 + a;(l-2;)fc2]5/2 ' 

The angular integrals are trivial and after introducing an 
ultraviolet cutoff A and rescaling the integration variable 
by I A: I the integral becomes 

C{k) = 27rfc2 / dx^ / dq^ ' ^—j^ . 

Jo Jo [q^ + x{l-x)]^'^ 

(B13) 

The remaining integrals are elementary. Isolating the 
leading infrared divergent term we obtain 



(B14) 



Substituting this to Eq. (|B9| ) we get the self energy (|6E|). 



3. Dimensionally regularized gauge field line 
integral 

To evaluate V{r) specified by Eq. we write it as a 
sum of two contributions, 



V 



8r 



2. TF self energy: Lorentz gauge 

For simplicity we evaluate the self energy (^) in the 
Lorentz gauge = 0). Extension to arbitrary covariant 
gauge is trivial. Eq. can be written as 



with 



^L{k)^-^J,I,ik) 



q-ik + q) 



I^,{k) ^ / d-'qq^ 



\q?{k + qY' 



(B9) 
(BIO) 



where we used an identity 



q,q. \ ^ _2 



Since the only 3- vector available is k, clearly the vector 
integral Ifj,{k) can only have components in the fc^ direc- 
tion, Ifj,{k) = C{k)k^/k'^ . By forming a scalar product 
kfilf_i{k) we obtain 

m . k,i,ik) = I d^ ^''\';l^':lP - (BID 



\q\Hk + qy 



where 



hir) 



„ik-r 



d^k- 



k ' 



ik-r ( L , ^2 



(B16) 
(B17) 



We first consider Ii{r). It is convenient to exponentiate 
the denominator by use of the formula 



1 



a-l-sA 



A'^ r(a) Jo 
Taking A = fc2, a = i we have 
1 r ds 



dss'^-'e 



h{r) = —= I —= I a- He 

Jo vs 



(B18) 



(B19) 



The k integral is now easy to evaluate by completing the 
square, shifting the integration variable, and making use 
of J d'^kexp{-sk^) = {ir/sy^^. We thus obtain 



d-l / d+1 _r£ 

/j(r) = TT 2 / dss 2 e 
Jo 



(B20) 
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Substitution t = /As tr ansfo rms (B20) to an integral 
of the type shown in Eq. (B18) and the result reads, 



Using the same procedure we obtain 

/2(r)=-(4^)^(d-2)r(^)ri-'^. 



(B21) 



(B22) 



Substituting ( |B2lD and ( |B22| ) back into ( piSD we obtain 
the result quoted in the text, Eq. (B3). 



APPENDIX C: PROOF OF BROWN'S RELATION 
FOR GAUGE INVARIANT PROPAGATOR 

To prove Brown's relation Eq. (|78| ) it is useful to in- 
troduce fermion propagator for a fixed configuration of 
gauge field a^. 



G[x, x'; a] 



1 



I?[T,T] T(x)T(x')e 



(CI) 



where 5'[T,T,a] = / (PrCo and Co is the QED3 La- 
grangian given by Eq. (|5|), z[a] = J V[T,T] e-'5[T,T,a] 
It is easy to see that G[x, x' \ a] solves the Dirac equation 



- a^)G[x,x';a] = 5{x - x'). 



(C2) 



Furthermore, in terms of G[ ] the TF propagator 

can be written as 



G{x-x') 



1 



Va G\x,x':a]e-'"' 



-SbW] 



(C3) 



where Sb = J (PxCb is the effective action [for our 
case, to one-loop £_b is given by Eq. (^] and Z = 
J Va e~'^«[''l . We now introduce a gauge invariant analog 
of G[x, x'; a], 



g[x,x';a] = e-'-f^' ''■''''G[x,x';a] 



(C4) 



where the integral in the exponent is taken along the 
straight line connecting space-time points x' and x. The 
gauge invariant TF propagator define d by Eq. ( [76| ) is then 
given by an expression analogous to (C3), 



g{x -^) = \ j'Da g[x, x'; ale-^^I"^! 



(C5) 



Our task is to relate Brown propagato r Q( x — x') to 
G{x~x') @. To this end we rewrite Eq. as follows 

G{x~x') = ^ J Vag[x,x';a\e'-'"-'^''^''^ 



4 / Va g[a;,a;';a]e-^^f''l 



g(x-x') 



(C6) 



where Snia] = Ssia] - iJ ■ a, Z = f Va e"-^^!"! and the 
last equality in Eq. (C6) should be taken as a definition 
of ^(x — x'). The source term J is defined by Eq. ( |80| ) 
and J • a is a shorthand for J (firJ{r) ■ a{r) = a ■ dr. 

Note that the linear UV divergence of phase factor in Z 
cancels out between numerator and denominator (|C6|) . 
We now observe that 

- = 1. I Va e''^"^e"^«M = (e*'^"^) = g-^^^"^'), (C7) 
Z Z J 



with F defined by Eq. 
we have (78): 



91). Substituting this into (|C6| ) 



G{x-x') = 



-F(a;-2;') 



Q{x-x'). 



(C8) 



To complete this part it remains to address the relation 
between Q{x — x') and Q{x — x'). To this end we no- 
tice that both Q and Q are gauge invariant: the former 
by construction and the latter by the following simple 
observation. Since both Q[x,x';a] and Ssla] depend on 
the transverse part of a, any dependence on longitudi- 
nal part of a comes through the iJ ■ a term and there- 
fore identically cancels between the numerator and Z in 
the denominator. Therefore, once defined by Eq. (C6) 
as a ratio of two gauge-variant objects computed in the 
same covariant gauge, Q takes on life of its own and is 
a fully gauge invariant quantity. We may rewrite G as 
such ratio in arbitrary linear gauge; specifically we may 
chose the axial gauge in which the iJ ■ a vanishes. In this 
gauge Q is exactly equal to Q. Since they are both gauge 
invariant they must be equal in arbitrary gauge. The 
reader should be cautioned, however, that this equality 
is a formal one since it involves manipulations of gauge- 
dependent quantities which are ill-defined in absence of 
some specific regularization. 



APPENDIX D: FEYNMAN INTEGRALS FOR 
ANISOTROPIC CASE 

1. Self energy 

As shown in the Section V, to first order in 1/N expan- 
sion, the topological fermion self energy can be reduced 
to 



Sn(g) 



d^k {q-kU2gl^i:-ilg-^;)Dx^,{k) 



(27r)3 



(DI) 



where Df^j, (q) is the screened gauge field propagator eval- 
uated in the Section 3. In order to perform the radial 
integral, we rescale the momenta as 



(D2) 
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and obtain 



(27r)3 



VFVAiQ-K)^ 



(D3) 



At low energies we can neglect the contribution from the 
bare field stiffness p in the gauge propagator (see Section 
3) and the resulting D^u{q) exhibits ^ scaling 



\K\ 



(D4) 



Now we can explicitly integrate over the magnitude of the 
rescaled momentum K by introducing an upper cut-off 
A and in the leading order we find that 



Repeated indices are not summed in the above expres- 
sion, unless explicitly indicated. In the case of weak 
anisotropy {vp = 1 + e, = 1) we can show that the 
Eq.(^) reduces to Eqs. (|ll5[]T]l ). 



Consider now the effect of the (covariant) gauge fix- 
ing term on rj^. Let us define the part of F^i, which 

depends on the gauge fixing parameter ^ as F^^K The 
general form of this term is Fj^f) = ^ kfik^f{k) where 
f{k) is a scalar function of all three components of fc^; 
f{k) does in general depend on the anisotropy. Upon 
rescaling with the nodal metric, see Eq.(D2), we have 



F, 



(?) _ 



sponding scalar function of K^. 
the Eq. (pSl) we find 



=Ki, f{K), where f{K) is the corre- 
Substituting Fj^^) into 



KiQ- 



Kf 



-'I 



Qx-2KxK-Q 
(D5) 



where K — {cos 9, sin 6* cos 0, sin 6* sin 0). Since is 
odd under inversion while F^^^, is even, it is not difficult 
to see that the term proportional to A vanishes upon the 
angular integration. Thus 



+Y.Ak^k,k^kX (D9) 



{2'kYvfva 



(Qx~2k^k-Q) [x] 



{x){2^,^Y:--lx9;.)F,Ae,cf,)\n (^^^ (D6) 

Using the fact that the diagonal elements of F^,y{9, (f) are 
even under parity, while the off-diagonal elements are 
odd under parity, the above expression can be further 
simplified to 



Sn(g) = -5](7;:9M^;:)in 



A 



(D7) 



where the coefficients rf^ are pure numbers depending on 
the bare anisotropy and are thereby reduced to a quadra- 
ture 



vfva{2tt)^ 



{2k, k^ - l){2g;^F,, - ^ 



-J24.k,k, 



^,,F^u]. (D8) 



where 77^ is the part of 77^ which comes entirely from the 

gauge fixing term. Using the fact that k^k^ = 1, it 
is a matter of simple algebra to show that 



4- 



dn f{K) 
VFVA{2Tr)^ 



(DIG) 



i.e., the dependence on the index ^ drops out! That 
means that the renormalization of 77^ due to the un- 
physical longitudinal modes is exactly the same for all 
of its components. Therefore, the difference in 771 and 7^2, 
which is related to the RG flow of the Dirac anisotropy 
comes entirely from the physical modes and is a gauge 
independent quantity. Note that this statement does not 
depend on the choice of the covariant gauge, i.e. on the 
exact form of the function /, only on the fact that the 
gauge is covariant. 
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